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rS ' Abstract 

We derive a factorization formula for the production of pairs of heavy coloured 
particles in hadronic collisions near the production threshold that establishes fac- 
torization of soft and Coulomb effects. This forms the basis for a combined re- 
summation of Coulomb and soft corrections, including the non-trivial interference 
of the two effects. We develop a resummation formalism valid at NNLL accuracy 
using the momentum-space approach to soft gluon resummation. We present nu- 
merical results for the NLL resummed squark-antisquark production cross section 
at the LHC and Tevatron, including also the contribution of squark-antisquark 
bound states below threshold. The total correction on top of the next-to-leading 
order approximation is found to be sizeable, and amounts to (4-20)% in the squark 
mass region 200 GeV - 3 TeV at the 14 TeV LHC. The scale dependence of the 
total cross section is also reduced. 



1 Introduction 

In perturbative calculations of partonic cross sections at hadron colliders there often arise 
terms that are kinematically enhanced in certain regions in phase space. In the Drell-Yan 
process, for instance, Sudakov logarithms of the form [a s In (1 — z)] n arise, where z = Q 2 /s, 
and Q 2 and s represent the invariant mass squared of the lepton pair and the partonic 
centre-of-mass (cms) energy, respectively. In the threshold region z — > 1, these terms are 
large, and spoil the convergence of the perturbative expansion in the QCD coupling a s . 
If the dominant contribution to the hadronic cross sections originates from the partonic 
threshold region, which is certainly the case when Q 2 approaches the cms energy s, these 
logarithms need to be resummed to all orders in perturbation theory to attain a reliable 
theoretical description. This was accomplished in [Tj[2] for the inclusive Drell-Yan cross 
section by solving evolution equations in Mellin space. In production processes of pairs 
of heavy coloured particles such as top quarks or coloured particles in extensions of the 
standard model, e.g. squarks and gluinos in supersymmetric extensions, the partonic cross 
section contains terms of the form (a s log (3 ) n ("threshold logarithms"), where — (1 — 
AM 2 /s) 1 ' 2 is the heavy particle velocity, and (a s /(3) n ("Coulomb singularities"), which are 
also enhanced near the partonic threshold s ~ AM 2 . Resummation of threshold logarithms 
for heavy-particle and di-jet production has also been implemented in Mellin space [3H9] 
up to now, and has been used for improved predictions of the top-pair production cross 
section at hadron colliders [10^16j. Recently the resummation of threshold logarithms 
for production processes of supersymmetric coloured particles [TTH2"U] and colour octet 
scalars [2"T||2"2"] has been studied as well. The resummation of multiple exchanges of Coulomb 
gluons and bound-state effects has been studied for the total top-antitop cross section [101 
[23] and the invariant mass distribution of top quarks and gluinos [21H26] . 

The theoretical basis for resummation is a factorization of the partonic cross section a 
in the partonic threshold region into hard and soft contributions of the schematic form 

a = H®S (1.1) 

with a hard function H and a soft function S both of which are matrices in colour space. 
If the invariant mass Q 2 = (pi + P2) 2 of the heavy-particle pair H(pi)H'(p 2 ) is held fixed, 
the "partonic threshold" is then defined more generally as the kinematical region where Q 2 
is close to the partonic centre-of-mass energy s. Arguments from perturbative QCD [27] , 
the properties of Wilson lines [2SI29], or effective theories [30H33] can then be used to 
demonstrate exponentiation of the enhanced contributions by solving evolution equations 
for the functions S and H. Traditionally this resummation is performed in Mellin space 
which requires a numerical transformation of the resummed result back to momentum 
space. A resummation method directly in momentum space was proposed in [21] and has 
been applied subsequently e.g. to Drell-Yan or Higgs production [33H35] and recently to 
the top-quark invariant mass distribution |36j. 

In this paper we consider the production of a pair of heavy coloured particles H, H' 
with masses m H and m H r, respectively, in the collision of hadrons Ni and N 2 , 

N 1 (P 1 )N 2 (P 2 ) -+ H( Pl )H'(p 2 ) + X (1.2) 



and concentrate on the resummation of the total cross section. In this case, after integration 
over the invariant mass of HH', the enhanced logarithms appear when 

8 « (m H + m w ) 2 . (1.3) 

In this limit, the factorization of soft gluons is complicated by the fact that the non- 
relativistic energy of the heavy particles is of the same order as the momenta of the soft 
gluons, in contrast to the assumptions made in the derivation of the usual factorization 
formula ( II. ip . We show that in this kinematical regime the partonic cross sections factorizes 
into three contributions, 

a = H(g)W(g)J, (1.4) 

where H is determined by hard fluctuations, W by soft fluctuations and J accounts for the 
propagation of the heavy-particle pair including Coulomb-gluon exchange. The short-hand 
W (g> J includes a convolution which accounts for the energy loss of the heavy particles due 
to soft emissions (the precise form of the factorization formula is given in (I2.63J) below). 
The main points of the present approach are as follows: 

• Our approach is largely model-independent and highlights the universal features of 
soft-gluon and Coulomb resummation: the precise nature of the heavy particles and 
the physics model enter the factorization formula (I1.4J) only in the hard function 
H; the soft function and the Coulomb function depend only on the colour and (in 
the latter case) the spin quantum numbers of the heavy-particle pair. The hard 
function can be obtained directly from the HH' production amplitude, expanded 
near threshold, without any need to perform a full cross section calculation. 

• Eq. ( II. 4p generalizes previous results of the form (II. ip by factorizing the Coulomb 
effects in addition to those from soft gluons. Multiple exchange of Coulomb gluons 
associated to corrections ~ a™//3 m can be resummed in the function J using methods 
familiar from non-relativistic QCD (NRQCD). The presence of the Coulomb function 
J leads to a more complicated colour structure of the soft functions W compared 
to previous treatments based on a factorization of the form ( II. ip . The factorization 
formula (jl.4j) allows a combined resummation of soft and Coulomb gluons and justifies 
earlier treatments [HKTHIEIIES] where the factorization of Coulomb from soft gluons 
was put in as an assumption. 

• The kinematical structure of the soft function W simplifies in the threshold re- 
gion (11. 3p compared to the more general situation considered in [TJ. This has al- 
lowed us to construct a basis in colour space that diagonalizes the soft functions 
relevant for hadron collider pair production processes of heavy particles in arbitrary 
representations of the colour gauge group to all orders of perturbation theory [37J. 
The diagonal colour bases for resummation of threshold logarithms for all produc- 
tion processes of pairs of coloured supersymmetric particles at hadron colliders, i.e. 
squark-antisquark, squark-squark, squark-gluino and gluino-gluino production have 



been provided in [37], extending explicit one- loop results for soft anomalous dimen- 
sions at threshold [7J[T7J[THJ[2n]- This result greatly simplifies threshold resumma- 
tions at NNLL accuracy and has allowed us to extract the two-loop soft anomalous 
dimension for arbitrary colour representations [37] from results of [38H4"0] . in agree- 
ment with an independent two-loop study for top-pair production [JT]. The more 
complicated colour structures in the two-loop soft anomalous dimension for massive 
particles with generic kinematics |4"2"H4"4] have been shown not to contribute to the 
NNLO total production cross section of a heavy-particle pair produced in an 5-wave 
at threshold 



We use the method of [31,32] to perform the resummation of threshold logarithms 
log/3 directly in momentum space solving evolution equations for the functions H 
and W . The resummation of the total cross section in the threshold region has been 
performed previously in Mellin space [9][TT]. Our derivation gives a field theoreti- 
cal definition of the quantities appearing in that approach and the resummation in 
momentum space allows an analytical treatment, since no numerical inverse Mellin 
transform is necessary. The relation between the Mellin-space and momentum-space 
formalisms has been discussed in 



The formalism includes the case of heavy particles with sizeable decay widths. To 
first approximation finite-width effects enter only in the Coulomb function and can 
be included by a shift of the non-relativistic energy E — > E + iT as familiar from 
top-quark pair production at electron-positron colliders. A systematic description 
of finite width effects beyond leading order can be achieved in the effective theory 
approach pf6TT49] . 



In this paper we focus on establishing the factorization of soft and Coulomb gluons and 
demonstrate the new features arising in a combined resummation by presenting numerical 
results for the example of squark-antisquark production. While we set up a method suitable 
for resummation at NNLL accuracy, in our initial application we restrict ourselves to 
NLL accuracy since the required colour separated hard production coefficients for squark- 
antisquark production are currently unknown. We also do not include finite squark decay 
widths which are small in typical supersymmetry (SUSY) scenarios and would introduce a 
dependence on the SUSY parameters but, as mentioned above, could easily be included in 
our approach. Note that, since in the present work we are concerned with the total cross 
section, there are no sizable finite-width corrections as long as V -C M, unlike the case of 
the invariant mass distribution near threshold. 

The paper is organized as follows. In section [2] we derive the factorization formula (I1.4J) 
using effective field theory methods. We compare to previous treatments of combined soft 
and Coulomb corrections and extend the diagrammatic argument of [H] on the absence of 
corrections from subleading couplings of soft gluons to non-relativistic particles at NNLL 
to all orders in the strong coupling using the effective theory approach. In section [3] we 
discuss the process-independent ingredients in the factorization formula, the soft function 
W and the Coulomb function J, collect their explicit results and perform the summation of 



threshold logarithms using evolution equations in momentum space. We also give a simple 
prescription for how to obtain the hard-function H, the only process-dependent ingredient 
in the resummation formalism, from a standard fixed-order calculation. In section H] we 
perform the combined soft-Coulomb resummation at NLL accuracy for squark-antisquark 
production and present numerical results for the Tevatron and the LHCo A number of 
additional results concerning parton densities, the Coulomb potential in various colour 
representations and the resummed and expanded cross section are collected in the appendix. 

2 Factorization for coloured heavy-particle pair pro- 
duction near the partonic threshold 

In this section we derive a factorization formula of the form ( 11. 4p for the process (11.20 near 
the partonic threshold. The heavy particles are assumed to transform under representations 
R and R' of the colour gauge group SU(3)c- The inclusive production cross section is 
described by the factorization formula 



a = 

p,p' 



^ / dXldx 2 fp/NtiXl, H)fp'/N 2 {x2, H) &pp'(xiX2S, fi), (2.1] 



where a pp ' are the hard-scattering cross sections for the partonic subprocesses 

p(ki) P '(k 2 ) -)■ H( Pl )H'(p 2 ) + X (2.2) 

and pp' G {qq,qq,qq,gg,gq,gq}- As usual the nucleon masses are neglected and parton 
momenta k\ t2 are related to the incoming nucleon momenta by 

h = Xl Pi , k 2 = x 2 P 2 . (2.3) 

The partonic cms energy is given by 

s = (h + k 2 f = Xl x 2 s. (2.4) 

At hadron colliders s is not a fixed quantity but one may argue that the steep rise of the 
parton distribution functions with decreasing x leads to an enhanced contribution to the 
total production cross section from the partonic threshold region 

AM 2 , s 

z = 1, (2.5) 

s 

M = (mu + m H>)/2 denoting the average heavy-particle mass. According to a recent 
study of threshold effects in the Drell-Yan process [33] this dynamical enhancement of the 
partonic threshold region remains effective for values 4M 2 /s = zxix 2 > 0.2. Therefore we 
expect the resummation of threshold logarithms to be relevant for particles with masses 



1 Preliminary results have appeared in [5U1I51J . 



M > 3 TeV (1.6,2.2 TeV) at the LHC at 14 TeV (7, 10 TeV) and M > 450 GeV at the 
Tevatron which are at the border of the estimated reach for coloured SUSY particles at the 
LHC [S^lol?] and the exclusion limits obtained at the Tevatron [54755] . However, experience 
with the NLO corrections to top-quark production suggests that also for smaller masses 
there is a sizable contribution to the total cross section from those terms in the partonic 
cross section that are enhanced for z — > 1 so that resummation can serve as an estimate of 
higher-order effects. In addition, expansions of resummed cross sections can also be used to 
predict the threshold-enhanced terms in higher-order perturbative calculations [T2~ } I41) } I5(3] . 
The heavy particles produced near threshold are non-relativistic. In terms of dimen- 
sionless variables the partonic threshold is characterized by the smallness of 



/ 4M 2 , . . 

/3 = yi- — = v / I^, (2.6) 

which for equal masses also corresponds to the non-relativistic velocity of the heavy parti- 
cles in the partonic cms. The aim of this paper is to perform a resummation of threshold 
logarithms log/3 and Coulomb corrections 1/(3. In order to discuss the systematics of the 
combined resummations of the two corrections we count both a s ln/3 and a s /(3 as quanti- 
ties of order one and introduce a parametric representation of the expansion of the cross 
section of the form |3j 



<7pp> OC <T 



' ^ ( -J ) ex P hi/3go(a 3 ln/^ +gi(a s ln/3) + a s g 2 (q s ln/^ +.. 
k= ° (LL) (NLL) (NNLL) 

x {1 (LL,NLL); a s , (3 (NNLL); a 2 s , a s /3 } /3 2 (NNNLL); . . .} . (2.7) 

With this counting, the resummed cross section at LL accuracy includes all terms of or- 
der l/(3 k x a™ +k In n (3 relative to the Born cross section near threshold. Next-to-leading 
summation includes in addition all terms of order a s ln(3; a 2 s {l/(3 x ln/3,ln 3 /3} ; . . ., while 
furthermore all terms a s ; a 2 {l//3, In ' /?};... are included in NNLL approximations 

In the threshold region, the cross section of the process (j2.2p receives contributions 
from hard, potential, soft and two collinear momentum regions defined by 



hard (h) 


k ~ M 




potential (p) 


k ~ MX, k ~ Ma/X 




soft (s) 


k ~ k ~ MX 




n-collinear (c) 


k_ ~ M, k + ~ MA, k ± 


~M\/A 


n- collinear (c) 


k + ~ M, k„ ~ MA, k ± 


~ M\/~X 



(2- 



2 The NNLL terms odd in f3 at order a s and a 2 sl a s /31n 2,1 /3; a 2 /31n 4:3 /3, vanish due to rotational 
invariance for the total cross section when the heavy-particle pair is dominantly produced in an iS-wave 
state at tree level [37N2S] . See also below in section 12.51 



assuming equal heavy-particle masses for simplicity. Here we introduced the generic power 
counting parameter A ~ 1 - z = /3 2 < 1. For the collinear momentum we employ the 
light-cone decomposition 

p^ = P- n ^ + P±n' 1 + P fl ± (2.9) 

with two light-like vectors n and n satisfying n • n — 2. Due to the threshold kinematics 
(k\ + &2) 2 ~ 4M 2 , no collinear modes can appear in the final state so that the particles 
denoted by X in (jl.2p are entirely given by soft modes. 

In order to achieve the separation of the effects corresponding to the different re- 
gions ( 12. 8 p it is useful to describe the heavy-pair production process in the framework 
of effective field theories. The relevant effective Lagrangian contains elements of non- 
relativistic QCD (NRQCD), describing the interaction of potential and soft modes, and 
soft-collinear effective theory (SCET) [5TH60] describing the interaction of soft and collinear 
modes. Hard modes with virtuality of order M 2 are integrated out and are not part of 
the effective Lagrangian. This effective theory is similar to that used in pi8lPi9l |6T] to 
describe W-pair production near threshold in e + e~ collisions and the construction of the 
effective Lagrangian proceeds along similar lines. As in the discussion of the Drell-Yan 
process in [33] the following considerations will be adequate to establish factorization of 
the partonic cross section; the discussion of the hadronic cross section requires to include 
additional modes with a scaling given by powers of Aq CD /M [32j. We also note the exis- 
tence of an additional region with scaling k ~ MvA, usually referred to as "soft" in the 
NRQCD literature. It is related to potential forces between the heavy particles, and will 
not be relevant for the following discussion of the leading Coulomb corrections. 

In the EFT framework, the factorization formula is derived by a series of steps involving 
first a matching calculation from QCD (or another "full theory" such as the MSSM) to 
the EFT followed by field redefinitions that decouple soft gluons from collinear 158062] and 
potential [63] fields. These will be described in the remainder of this sectionjj We also 
show that, up to NNLL accuracy, corrections from sub leading interactions to the effective 
Lagrangians not included in the derivation of the factorization formula either vanish or can 
be straightforwardly incorporated. 

2.1 Matching to SCET and NRQCD 

The first step in the derivation of the factorization formula consists of representing the 
scattering amplitude for the process (11.21) in terms of effective theory matrix elements of 
production operators C{ a . a i corresponding to the partonic sub-processes p(ki)p'(k 2 ) — > 
H(pi)H ; (p2)X, multiplied by coefficient functions that contain the hard part of the ampli- 
tude. After this step the partonic scattering amplitude is given by 

A(pp' -> HH'X) = J2 CfaU^ (HH'X\Ofl a} (ri\pp') EFT . (2.10) 



3 Thc corresponding discussion of factorization of electromagnetic effects and resummation for the W- 
pair production cross section at e + e~ colliders has been given in |64j . 
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The subscript "EFT" means that the matrix element is evaluated with the effective La- 
grangian discussed below. We now explain the definition of the operators and the calcula- 
tion of the coefficient functions. 

The hard production subprocess corresponds to the 2 — > 2 process p(ki)p'{k,2) — > 
H{pi)H'{j)2)- Near the partonic threshold, kinematics forbids additional massive or en- 
ergetic particles in the production operator, while additional soft fields connected directly 
to the hard process imply additional highly off-shell propagators and therefore lead to op- 
erators power-suppressed in A ~ 1 — z relative to the leading four-particle operators. The 
required four-particle operators take the form 

<1}M = [^a^-o^AsaA,] (A*)- (2-11) 

The form of the production operators is thus similar to those describing the production 
of non-relativistic H^-pairs in e + e~ collisions discussed in |18|H9]. Several remarks on the 
notation are in order: 

(1) The fields ip^ and ip'^ are non-relativistic fields which create the heavy particles H 
and H' with momenta 

Vi — rfiuw^ + p[ , T>2 = m fl' w '' + P2 (2-12) 

with the cms velocity w^ (w 2 = 1) and small residual momenta pi/2 in the potential region 
(12. 8p . In the partonic cms frame w = (1,0) and p% = — p~2 = P for the spatial components 
of pi^- The fields <p c (0 g ) are collinear (anticollinear) fields that destroy the initial state 
partons p and p' with momenta 

fef « Mn» , fc£ « Mff, (2.13) 

respectivelyjj We provide more details on the definition of these fields and their effective 
Lagrangians in the following subsection. 

(2) Greek indices oij denote the spin (or Lorentz group representation) index of the 
field. Where convenient we use the multi-index notation {a} = a\a2a^a\. Repeated 
multi-indices are summed by summing over all four components. A similar convention 
is used for colour indices, which we denote by Latin letters {a} = a,\ . . . a 4 . We employ 
creation fields for the heavy final state particles and destruction fields for the initial state. 
If, for example, H' is the antiparticle of H, then this convention implies that ip' transforms 
in the complex conjugate SU(3) representation of ip. Similarly, an initial-state antiquark 
transforms in the antifundamental representation. In effective-theory calculations it is 
conventional to decompose operators into a complete basis in spin and colour, making 
operators and coefficient functions "scalars". Here we prefer to work with operators and 
coefficient functions that carry open spin and colour indices, see (12. lip . This turns out to 



4 We do not have to consider operators where the (f> c destroy particles different from p and p' . These 
would arise from a splitting of a collinear gluon into a collinear and a soft quark or of a collinear quark 
into a soft quark and a collinear gluon, but these interactions are present only in the power-suppressed 
part of the SCET Lagrangian [59] . 



be convenient for the spin indices, since the soft gluon contributions are spin-independent, 
as well as the Coulomb contributions needed at NLL accuracy. Thus, it is more direct to 
perform the spin summations on the square of the amplitude as in standard unpolarized 
cross section calculations instead of decomposing the amplitude. At NNLL accuracy the 
spin-dependence of the potential corrections has to be taken into account in order to 
obtain all terms of the order a 2 s In /3. For this purpose a spin decomposition of the potential 
function is introduced in section [231 As concerns colour, given the representations r, r' and 
R, R' of the initial and final state particles, we can introduce a set of independent colour 
matrices c- 1 ' that form an orthonormal basis of invariant tensors in the representation 
r £g> r' <g> i? <g> R ' . The colour structure of the Wilson coefficients can be decomposed in this 
basis according to 

i 

This allows us to write 

c£.,m o£«>m = E c t>) [4Wt»w] - E cg?<*o °l>) (2-15) 

i i 

in (12.101) . For example, in a 3 + 3 — > 3 + 3 scattering process (such as qq — » ti and 
squark-antisquark production), a convenient colour basis is [7] 




c (1) - 

(2) _ z T b T b 
{a} /TT~ «2"1 tt 3»4 ' 



(2.16) 



with Da = N% — 1, expressing colour conservation and the fact that the initial and final 
states can be in a colour- singlet or a colour-octet state. Suitable colour bases for the 
different colour scattering processes have been constructed in [37J where also more details 
on the properties of the basis tensors can be found. 

(3) In (J2.10P the argument // denotes the factorization scale dependence of the operators 
and hard functions. We drop this in the following, but note that exploiting the invariance 
of the physical amplitude under changes of the factorization scale will be a key element 
in deriving the resummation formula. The superscript C. stands for higher order terms in 
the expansion of the amplitude, which correspond to operators of the form (12. lip with 
derivatives acting on the fields, as discussed below. 

(4) Fields without space-time arguments are at x — 0. In general, the SCET repre- 
sentation of a scattering amplitude involves a convolution of the Wilson coefficient with 
collinear fields evaluated at different positions along the light ray [59], instead of a prod- 
uct with a local bilinear operator 4> C ;a 1 a 1 ( t>c;a2a2 as i n (I2.10p . However, since here there 
is only one collinear field of a given type in the operator, translation invariance may be 
used to rewrite the operator in a local form [591165] . Applying the translation operators 
on the initial state turns the convolution integrals with the coefficient function into the 
Fourier-transform of the position-space coefficient functions to momentum-space coefficient 
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functions. These then depend on the large momentum components of the partons, and 

y {a;a} 



are the natural objects to work with. The hard functions CP, in (I2.10p refer to the 



momentum-space quantities. 

Despite the somewhat complicated formalism and notation, there is a straightforward 
general prescription for the computation of the matching coefficients C, , in the expres- 
sion ( 12.1 Op . One simply calculates the renormalized scattering amplitude for the par- 
tonic process p(ki)p'(k2) — > H(pi)H'(p 2 ) without any averages or sums over colour and 
spin states, and expands it in p^ 2 defined by (12.1 2[) around the production threshold 
ik\ + k 2 ) 2 = (mfl + niH') 2 - Beyond tree level, the matching coefficients can be extracted 
directly by expanding the loop integrand under the assumption that all loop momenta are 
hard [66J and by performing the loop integration of the expanded integrand in dimensional 
regularization. As a special case, the leading term in the kinematical expansion is found 
by calculating the full amplitude directly at threshold. In this case the loop corrections to 
the effective theory matrix element evaluated directly at threshold vanish since they are 
given by scaleless integrals. Therefore only tree diagrams have to be considered on the 
effective theory side of the matching calculation and the hard function is determined from 

A(pp' -> HH% =4KP = J2 Cfi} (HH'\0^. a} \pp'} EFTtree . (2.17) 

i 

Since the matrix element on the right-hand side is simply a product of spinors (polarization 
vectors) and the operator renormalization factor, the MS scheme coefficient function CJJ a -i 
equals the scattering amplitude at threshold with the spinors (polarization vectors) stripped 
off, and all 1/e poles set to zero. The Born- level matching required for squark-antisquark 
production is discussed explicitly in section HI 

The leading term in the kinematical expansion describes the production of HH' in an 
S- wave state. The next term is related to the terms linear in p 12 , more precisely their 
components orthogonal to u> M , denoted by the symbol T, since the other component is of 
higher order in the potential region. The P-wave production operators corresponding to 
this term are of the form (12.1 ip with the replacement 

(2.18) 
in the non-relativistic part of the operator, with the covariant derivative D^ in the appro- 
priate colour representation. The series in C. in ( I2.10p accounts for the threshold expansion 
in powers of A. In the following, we will mostly consider only the leading term £ = 0. In 
particular, the factorization formula we derive is valid only for the case that the Born cross 
section is dominated by S'-wave production. 



2.2 Effective Lagrangians 

2.2.1 Non-relativistic 

The propagation and interactions of the non-relativistic fields ip and if}' are described by the 
potential NRQCD Lagrangian [67H7U] . That is, we assume that the modes with momentum 
k ~ MyX and potential gluons have already been integrated out. At leading order, in 
the partonic cms frame, the effective Lagrangian, including the decay widths of the heavy 
particles, is given by 



C 



PNR.QCD 



Vit[iD? + _?- + »^U + vi'tf<i?» + ^- + ,T "' 



d 3 f 



2m H "2 
^T {R)a ip 



'a. 



2m h' 

V/t T («')<y] (0) . 



i// 



(2.19) 



Here T^) are the SU(3) generators in the appropriate representation of the heavy par- 
ticle. In this convention if>' transforms in the antifundamental representation, if H' is an 
antiparticle in the fundamental representation. This deviates from the standard NRQCD 
notation for quark-antiquark production that employs a field \ that creates an antiquark 
and transforms in the fundamental representation. Our notation is convenient since it 
accommodates the case where the heavy particles produced are not particle-antiparticle 
pairs. This convention leads to a different sign of the potential term compared to the 
conventional formulation. 

The only relevant interactions in ( I2.19P are the exchange of Coulomb (potential) glu- 
ons, rewritten as an instantaneous, but spatially non-local operator in the second line, and 
an interaction with the zero component (that is, w ■ A s , in a general frame) of the soft 
gluon field in the soft covariant derivative iD®ip(x) = (id + g s T^ a A a0 (xo))ifj(x), and the 
analogous expression for ip' with the generator in the representation R' . Using the canoni- 
cally normalized NRQCD Lagrangian ( I2.19p . one has to take into account a normalization 
factor V2E^ = {2m H {l + p 2 /m 2 H )) l l 2 = ^2m^ + 0(X 2 ) in the definition of the state \H) 
(similarly for \H')) on the NRQCD side of the matching relationjj 

The NRQCD Lagrangian (12. 19[) is invariant under soft gauge transformations in the 
appropriate representation: 



^{x) ->uW( Xq )i/,( x ) 



^{x)^Uf\x )^'{x). 



(2.20) 



The soft gauge transformation depends only on x since soft fields have to be multipole 
expanded when multiplied with potential fields in order to maintain a uniform power 
counting [70TT72"] . The invariance of the Coulomb potential for arbitrary representations 
follows from the property 

y(B)t T (fl)ay(fl) _ y(8) T (fl)6 (2.21) 



5 For Th — 0. See [47] for the general case. 
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where U^ R ' is a SU(3) transformation in the representation R and "8" denotes the adjoint. 
Therefore under a soft gauge transformation: 

[^t T («)<Y] ( x + r -) ^ fyt T (*)fy] ( x + F )U^ b (x ) (2.22) 

and analogously for the ip' fields. The gauge invariance of the Coulomb potential term 
follows then from the fact that the transformations in the adjoint representation are real 
so that U ab Uac = (U^U^)bc = Sbc- For this is essential that only the time component 
enters the soft gauge transformation since the ip fields and the ip' fields in the Coulomb 
potential are defined at different points in space but at the same time. 

2.2.2 Soft-collinear 

The propagation and interactions of quark and gluon collinear modes with large momentum 
proportional to n M is described by the SCET Lagrangian [ST]. At leading order, in the 
position-space formalism of SCET [59j[60] it is given by 

£c = & (in-D + % jp^-± w ^^ | £ c - i t r {F?F%) . (2.23) 

Here £ c denotes the n-collinear quark field, which satisfies r/£ c = and the projection 
identity (r/?^/4)£ c = £ c . The covariant derivatives are defined as iD c = id + g s A c with A c 
the matrix-valued gluon field in the fundamental representation. The covariant derivative 
D without subscript, however, contains both the collinear and soft gluon field. The quantity 
F^ u is the field-strength tensor built from the collinear gauge field in the usual way, except 
that n_ ■ D rather than n_ ■ D c appears [60]. The collinear and soft fields are evaluated 
at x, but in products of soft and collinear fields the soft fields are evaluated at x+ = 
(n ■ x/2) n M = x_n M , according to the multipole expansion. In this notation X- is a scalar, 
while x+ is a vector. The fi-collinear fields are described by an identical Lagrangian with 
the roles of n and n interchanged. The corresponding quark field satisfies ?^£ 5 = and 
(?^/4)£ g = £ a . The two collinear sectors couple only via soft gluon interactions. 

It is convenient to express the collinear part of the operators (12. lip in terms of fields 
that are invariant under the collinear gauge transformation as defined in [HD] ■ Let W c be 
the collinear Wilson line in the n direction, 



W c (x) = Pexp 



o 



ig s \ dtn ■ A c (x + nt) 



(2.24) 



Then (p c in (12. lip is given by the combinations 

W&{x), LW c (x) (2.25) 

for the quark and antiquark initial state, respectively. Note that our conventions imply 
that the conjugate fields (p\ are £} C W C = Ccl°W c for the quark and W^°(, c for the antiquark 
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initial state. For the gluon initial state we use the following definition that is invariant 
under the collinear gauge transformations 



r° 

Aj"{x) = g-\Wl[iD^W c }){x) = dsn u (W}F^W C ) (x + sn) . (2.26) 

J — oo 

The anticollinear fields are defined similarly with n replaced by n. These collinear operators 
transform under the soft gauge symmetry of the SCET Lagrangian [60] according to 



(W&)(x) -> UJ^(xJ){WiQ{x) 
A a c (x)T a -> Uf\x_)A a c {x)T a U^\x-) = U^{x.)A h c {x)T a 



(8), . .,. ( 2 - 27 ) 



t(R) 



where Us i.xJ) is a soft gauge transformation in the representation R at the point x^n^. 

2.3 Decoupling of soft gluons 

Since the collinear fields C , anticollinear fields g and the potential fields ip^ in the produc- 
tion operators (12. lip interact with each other only via exchange of soft gluons, an essential 
step in deriving a factorization of the hadronic cross section is the decoupling of soft glu- 
ons from the collinear and potential degrees of freedom. In SCET it is well known that 
the soft gluons can be decoupled from the collinear fields at leading power by performing 
field redefinitions involving Wilson lines [58] . Here we will show that an analogous trans- 
formation also decouples the soft gluons from the potential fields at leading order in the 
non-relativistic expansion. As a result, the production operators factorize into products of 
non-interacting collinear, anticollinear, potential and soft contributions. 

In order to decouple the soft gluons from the collinear and anticollinear fields, the 
required redefinitions are familiar from the derivation of factorization formulas in SCET [581 



ux) = g°^)si 8)t (s-) = ao^pv) 

K>"(x) = 0^)4* ( 2 - 28 ) 

with the soft Wilson lines Sn for a particle in the representation R of SU(3) given by 

r-0 



Sf) (*) =Pexp 



>fj: 



J dtn-A c s (x + nt)T {R)c 



(2.29) 



and similarly for particles collinear to n. Consistent with our treatment of incoming an- 
tiparticles as particles in the complex conjugate representations, we have rewritten the 
transformation of the conjugate collinear spinor as a transformation in the antifundamen- 
tal representation with the generators T^ = —T T . 
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These transformations induce analogous transformations on the collinear fields <p c G 
{Wj£ c , £ C W C , A}:} entering the production operators: 

<Pc-, aa (x) = S&ix-rfXix) (2.30) 

as follows from the transformation 

W c (x) -► S®{x-)w2>\x)S<M{x-) (2.31) 

of the collinear Wilson line. For the incoming parton in the anticollinear direction we 
exchange n by n, x_ by x + in the argument of the Wilson lines, and the representation 
r by r'. In the expressions with the superscript (0) all collinear fields are replaced by the 
decoupled fields A c and £c and the soft gluon fields are set to zero. 

Turning to the non-relativistic sector, we first note that from the diagrammatic per- 
spective the coupling between the non-relativistic and soft fields is non-trivial, since the 
non-relativistic energy of the heavy fields and the energy of soft gluons is of the same order 
MX ~ M/3 2 . Thus, soft gluons can attach to the non-relativistic lines, to and in between 
the Coulomb ladder rungs, leaving potential lines in the potential region. However, a redef- 
inition of the non-relativistic fields by a time-like Wilson line decouples the soft gluons also 
from the potential fields at leading order in the non-relativistic expansion. Namely, the 
interaction with soft gluons in the PNRQCD Lagrangian given in (I2.19P can be eliminated 
through 

Ux) = S ( w R) ab (x )4°\x) (2.32) 

ri(x)=il;W(x)S™(x ), (2.33) 



where the soft Wilson lines Sw are defined a; 
S { w R \x) = Pexp 
S^(x) = Pex P ,<j. 



oo 



/■oo 

/ dtw-A c s {x + wt)T {R)c 
Jo 

POO 

/ dtw A c s (x + wt) T {R)c 
Jo 



(2.34) 
(2.35) 



and analogously for the primed field transforming under the representation E! . This trans- 
formation eliminates the interaction contained in the soft covariant derivative D° s in (12.191) . 
since 

SWiBPSW = id , (2.36) 



6 Note that in contrast to (J2.30I) we do not use a field redefinition of the non-relativistic annihilation 
fields ip with Wilson lines extending from -co to ieo but (|2.32p . since the former definition would obscure 
the decoupling of the soft gluons in the PNRQCD Lagrangian. As discussed in [74], the two different forms 
of the redefinition are equivalent if one takes an appropriate phase for external non-relativistic in-states 
into account. Since we encounter only external outgoing non-relativistic particles, this subtlety is irrelevant 
in the present context. 
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but at first sight introduces the soft Wilson lines into the Coulomb potential interaction. 
However, since the identity (I2.2ip also holds if the group elements U are the Wilson lines 
S w [31], and since the S w in the transformations (I2.32p . (I2.33P only depends on the time 
coordinate, the soft gluon Wilson lines drop out from the Coulomb potential term ex- 
pressed in terms of the redefined fields. In other words, in terms of the redefined fields the 
PNRQCD Lagrangian takes exactly the same form as f 1 2 . 1 9 j) except that D° s — > d°, so that 
the Lagrangian is independent of the soft-gluon field. At the next order in the velocity 
expansion, the PNRQCD Lagrangian contains an interaction involving the gauge invariant 
chromoelectric operator x ■ E s (t,0) [67J[75] that is not removed by the above transforma- 
tion so the decoupling of soft and potential modes is only valid at leading order in the 
non-relativistic expansion!] We demonstrate in section [231 that the corrections to the total 
cross section from a single insertion of this operator vanish to all orders. 

Assembling the different contributions shows that the operators (12. lip for 5-wave pro- 
duction (£ = 0) factorize into an expression of the form 



)(0,0 („ „\ _ cd) 



with the universal soft contribution 

S{j}{ X ) = C {a}^n,a 1 j 1 {X-)Sn,a2J2{ X +)^w,j 3 a 3 { X 0)^w,j 4 a4,\ X 0) ■ (2.38) 

Here and in the following we suppress the representation labels on the soft Wilson lines, and 
take it as understood that the Wilson lines always are in the representations appropriate for 
the respective particles. The decoupled fields ijj^\ <fi c and (f) 5 in (I2.37J) do not interact 
with the soft Wilson lines and with each other at leading power in SCET and at leading 
order in the non-relativistic expansion. 

Near the partonic threshold only soft radiation can occur in the final state so the state 
\X) is free of collinear radiation. Therefore the Fock space is a direct product of potential, 
soft, collinear and anticollinear contributions, 

\PP) = \V) C ® \v')- c ® |0) s ® |0) p , \HH'X) = |0) c ® |0) 5 ® \X) S ® \HH') p , (2.39) 
and we obtain the scattering amplitude in the factorized form: 

A(pp'^HH'X)=J2 CgjV) (^1^(0)10) (M'lCVSJlO) 

i 

x(O|0S lQl b)(O|0g 2Q2 b')- (2-40) 

As remarked above, soft-potential and soft-collinear factorization is only valid at leading 
order in the non-relativistic and collinear expansion. Similarly for a P-wave production 
operator (I2.18P the field redefinitions (12.331) do not eliminate the soft gluons in the covariant 
derivative. 



7 In order to achieve factorization of soft gluons in a kinematical regime where relativistic corrections 
become important one can describe the heavy fields by two copies of heavy quark effective theory with 
different velocities wi and w%, see also [36]. Performing separate field redefinitions for H and H' with 
Wilson lines defined in terms of the two velocities decouples the soft gluons up to corrections of the order 
l/niH, I/tuh'- 
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(Tpp'{S) 



2.4 Factorization formula for the cross section near threshold 

We now show that the total partonic cross section for the production of a heavy-particle 
pair near the production threshold factorizes into hard, potential and soft contributions 
as sketched in ( II. 4p . To this end, we consider the cross section for the production of a 
heavy-particle pair near threshold in an S'-wave from the partonic process (12. 2p for on- 
shell, massless initial state partons p and p' with zero transverse momentum. In this case 
only the hard, soft, collinear and anticollinear modes (12.81) contribute to the cross section 
since no smaller scale related to Aqcd, parton masses or off-shellness is introduced and 
the modes included in the effective theory introduced in section [2TT1 are sufficient. We can 
therefore use the expression for the amplitude (12.401) in terms of decoupled collinear and 
potential fields and soft Wilson lines to compute the cross section as usual by squaring the 
scattering amplitude, averaging over initial state and summing over final state polarizations 
and colours, and integrating over the final-state phase space: 

j- J d^d^ 2 d^ x J2Yl \APP'^HH'X)\ 2 (2n) A 5^\p 1 +p 2 + P x -P) (2.41) 

pol colour 

Here P is the total incoming momentum and Px is the total momentum of the hadronic 
state X. We suppress the normalization factors for the initial-state averages. We do not aim 
here at a proof of the factorization of the hadronic cross section, that would also require the 
treatment of modes related to the QCD scale, as in the case of deep inelastic scattering for 
x — > 1 in [32] . Rather, we are concerned with the dimensionally regularized partonic cross 
section (12.411) and its short-distance counterpart with the initial-state infrared divergences 
subtracted. 

The factorized form (12.401) of the scattering amplitude holds if the expectation value of 
the EFT operator is evaluated with the leading order PNRQCD and SCET Lagrangians. 
This is a priori appropriate for NLL accuracy in the counting (12. 7p . Since according to 
(12. 7p the NNLL approximation includes 0(f3) corrections, subleading soft interactions in 
the effective Lagrangian that are not removed by the decoupling transformations could 
contribute to the cross section starting from NNLL. At fixed-order NNLO accuracy it was 
shown in |45] that these corrections do not, in fact, contribute to the total cross section and 
the only corrections not generated by the leading order effective Lagrangian are the one- 
loop hard function and NLO Coulomb and non-Coulomb potentials. These considerations 
will be elaborated on in section 12.51 and extended to all orders in the coupling expansion. 
Therefore the derivation of the simple factorization of the form (II .4p given in the following 
holds at NNLL accuracy for production processes dominated by S'-wave production. 

As for the hadronic cross section (12. ip . standard QCD factorization implies that the 
cross section for the parton initial states factorizes according to 

Vpp'(s) = ^2 dxidx 2 f P /p(xi, p)fp>/p'(x 2 , p) d- P p'(x 1 x 2 s, //), (2.42) 

p,p' 

where the parton distribution functions of parton p in parton p, f p /p(x,fi), contain the 
collinear modes. We will consider this cross section near the production threshold, s ps 
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4M 2 , and show that the short- distance cross section appi itself factorizes into hard, soft and 
potential contributions as sketched in ( II. 4p . At threshold, the PDFs in (j2.42p necessarily 
appear in the limit x — > 1, where the flavour off-diagonal PDFs are suppressed, so it is 
sufficient to consider the case p = p, p' = p' . This is consistent with the fact that, at 
leading power, we only have to consider EFT matrix elements in (I2.10p where the collinear 
(anticollinear) field C (0 g ) annihilates the state \p) i\p'))- According to QCD factorization, 
appi is independent of the initial states p, p' so it is identical to the hard-scattering cross 
section to be used in the hadronic cross section (12.11) . Therefore the factorized form of 
the short- distance cross section (ll.4p derived in this way can be used in the hadronic 
cross section in the partonic threshold region where X1X2S ~ 4M 2 . As discussed above, 
this region is expected to give a sizable contribution to the total cross section for heavy- 
particle masses 4M 2 ~ 0.2s. For values of X\p, outside the partonic threshold region, 
the factorized formula is not necessarily a good approximation to the cross section and 
the resummed cross section should be matched to the fixed-order expansion as discussed 
further in section [3T6l 

We now proceed with the derivation of the factorization formula which follows similar 
steps as that for the Drell-Yan process in [33]. The new ingredients here are the presence 
of the potential fields and additional Wilson lines for the heavy coloured particles, result- 
ing in a more involved colour and spin structure. After inserting the factorized matrix 
element (12.401) into the standard expression for the cross section (I2.4ip . we use an integral 
representation of the 5 function, 



(27r)^( Pl +p 2 + P x -P) = J d 



>4 ze -iz< P i+P2+Px-P)^ (2.43) 



Rewriting the exponential in terms of the momentum operator acting on the external 
states, we can translate the collinear, potential and soft operators in the conjugate matrix 
element to z, e.g. e iz ' kl (p(£;i)|4 0)t (0)|0) = (p(A;i)|e^4 0) *(0)e- iz ' # |0) = (pih)^ (z)\0) . 
We then obtain for the cross section: 

„ ( a \ — ^ V^ V^ V^ r (0,i) r (0,i')* I jA vf> -2iMwz jfafaazaif \ 

a pp'\ s ) - y s 2-^ 2_> 2^ °{«} u m J aze J k 3 k 4 j 3ji y z > 

pol colour i^i 

x WSU W|0> WgU(z)|0> (O|0g iai (O)b) (O|0g 2Q2 (O)b'> 

x jd^ x (0|<S { ( ;]»|X) (X\S® } (0)|0> . (2.44) 

Here we introduced the definition of the potential function in terms of the PNRQCD matrix 
elements: 

pol, colour 



Jffiiz) =<&«&* J2 (0|[^^£ ) a 1 ]Wl^fl'><^ff / |[^£a i ^22I](0)|0>. (2.45) 



The factor e 2tMw ~ z = e l( , m H+m H ,)w-z ar j ses because, by definition, the non-relativistic fields 
depend only on the residual momentum so that ip{z) = e zz ' p ip{0)e~ iZ ' P = e~ tz '^ip(0). We 
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can introduce a colour basis for the potential function J by introducing projectors P Ra 
onto the irreducible representations R a appearing in the decomposition of the final-state 
representations, R® R! = J2r R&'- 

J M=E P M J £ } - (2-46) 

Ra 

A systematic procedure for the construction of the projectors for arbitrary representations 
out of Clebsch-Gordan coefficients for the coupling R® R' — >■ R a is reviewed in [37] where 
the projectors needed for all production processes of pairs of coloured SUSY particles are 
collected and further properties of the projectors can be found. 

Since the leading order PNRQCD Lagrangian is spin independent, the spin structure 
of the leading order potential function (resumming multiple Coulomb exchange) is trivial: 

J £ } ( x ) = <W<5<*2a4 JrM- (2.47) 

However, as mentioned above, spin-dependent non-Coulomb potentials become relevant 
starting from NNLL accuracy in the combined counting in a s In /3 and a s /(3 defined in (12.71) . 
Therefore we introduce a decomposition of the potential function into different spin states: 

4*>)= E n ? } ^»- (2-48) 

S=\s-s'\ 

Here s (s') is the spin of the heavy particle P (P'). The projectors Us project on the HH' 
state with total spin 5*. They satisfy the completeness relation 











/ _j \a\ *aiazVa20L4 ■ 


(2.49) 










S=\s-s'\ 




As 


an example, 


for two 


spin 


| particles, the projectors on singlet and triplet states read in 


the 


cms frame: 






rrW - Ia A 

■ L1 r . u ct2aiU a3a4: 

rrM - 1^ cj'i- 

1 n a2«i 0304' 


(2.50) 



The explicit expression for the function J is given in section | 

In order to bring (I2.44p into the form of the standard factorization formula (12.421) we 
have to identify the collinear matrix elements with the PDFs which requires to combine the 
products of the matrix elements of the collinear fields into a single expectation value (the 
same discussion applies for the anticollinear fields) . This can be achieved by recalling that 
due to the threshold kinematics collinear radiation from the initial state is inhibited. We 
can therefore formally sum over the states of the collinear final state Fock space \C) since 
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the only contribution comes from the vacuum. After subsequent use of the completeness 
relation we can simplify the product of collinear matrix elements to 

(p\^(z)\0) (O|0f (0)|p) =► J2 [ d ®c (P\^^)\C) (C|0i o) (O)b) 

c J 

= (p|^ 0)t Wi 0) (0)b>- (2-51) 

The same result for the collinear matrix element is obtained in the standard collinear 
factorization away from threshold where the final state contains collinear radiation and the 
sum over the collinear Fock space is present from the beginning. We therefore define 

Jo X\ 

(2.52) 

where p is the parton annihilated by the field <fi and the average refers to colour and 
polarization. We show in appendix |A] that the so defined functions f p /p{xi, //) coincide with 
the standard definition of the quark, antiquark and gluon parton distribution functions. 
Since the momentum of the external state is Pf = n • P\ n M /2, the matrix element depends 
on z only through z + = n ■ z oc z ■ P\ since z ■ Pi = (n ■ z)(n ■ P\)/2, which has been used 
to represent it as a one-dimensional Fourier transform. This also allows us to replace z M 
by (n • z) W 1 /2 on the left-hand side. For all parton species, the prefactor is given by the 
helicity average of the polarization wave functions Xa of the SCET fields 

where the sum extends over the physical polarization states. D r denotes the dimension of 
the SU(3) colour representation of the initial state parton p and the number of spin-degrees 
of freedom is n s = 2 for all relevant cases. Note that the polarization wave functions are 
evaluated for a particle with momentum X\P\, which corresponds to the standard rule for 
calculating parton scattering cross sections. At the same time the factors l/x\ (and l/x^ 
for the other parton) from the definition (12.521) combine with the prefactor l/(2s) in (I2.44p 
to yield the standard flux factor l/(2s) for partonic scattering. 

The soft matrix elements in the cross section (12.441) can be collected in the soft function 

x J 

which is a matrix in the colour basis given by the c^K. Since the collinear matrix elements 
are diagonal in colour and due to the colour decomposition of the potential function (12.461) . 
the initial-state colour indices ji, ji are summed over and the final state indices are pro- 
jected on the irreducible representations R a . Since we are concerned with the total cross 
section, the soft matrix elements can be combined in a correlation function by summing 
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over the hadronic final state and using the completeness relation 

J2 fd$ x \X) (X\ = I (2.55) 

x J 

of the soft Hilbert space. The soft function can then be written as 

W«"(z, n) = P { f } cf a} W { { k J } (z, „)($> (2.56) 

with the correlation function of soft Wilson lines 

In combining the soft matrix elements into a correlation function we have introduced 
time- and anti-time-ordering symbols as discussed in [33J. In evaluating this correlation 
function, a soft gluon propagator connecting fields in the time-ordered and anti-time- 
ordered products is given by a cut propagator, therefore this prescription reproduces the 
usual rules for the real soft corrections. A colour basis that diagonalizes the soft matrices 
W w a to all orders in perturbation theory has been constructed in [37] as reviewed in 
section 13.11 below. 

Inserting these definitions, the cross section becomes a convolution of the PDFs with 
the potential and soft functions: 



Opp' 



dx 1 dx 2 fp/p(x ll (i)f I //pr(x2,lJt) Yl Yl H ii'^) 

i,V S=\s-s'\ 

x [w^i J L(l) I ' #zJ( Xl *-* a *>- 2Mi °-*)- z W$r(z,ti), (2.58) 

where we have introduced the Fourier transform of the potential function 

J S rM = J^^JrM ( 2 - 59 ) 

and the spin and colour dependent hard functions defined in terms of the short- distance 
coefficients averaged over the initial state spins and projected on the spin state of the 
heavy-particle system: 



1 



2.s 



(°>*) t ,,MT/ 3 3/34a 3 a4 / ^(0,i')* , 



HM = ~ K^Pi)K,^ p i) c Z'^r m "c y Z "(p.) 



= him: E^MxixSjiiS*-"^^)^^)' ■ <2 ' 60) 

A1A2 

In the second line we have inserted the prefactors (12.531) (omitting the momentum argu- 
ment of the polarization wave functions) Jj Summing over the spin states and using the 



8 In (|2.60p we could replace s by 4M 2 to make the hard function energy-independent. However, in the 
numerical evaluation it is trivial to keep the kinematic factor l/(2s) and this is the convention we adopt 
in Section [U 
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completeness relation of the spin projectors we obtain the spin-averaged hard function that 
is sufficient for NLL resummations: 

s-\-s' 

HM= £ ^(A*)=^[^ A (x 1 A)^(^2)C^(A*)^g&0')] (2-61) 

S=\s-s'\ 

A simple prescription for the computation of the hard functions directly from the scattering 
amplitudes will be given in section 13.21 

Comparing ( I2.58P to the factorization formula ( I2.42p we obtain the factorized expression 
for the short-distance cross section by stripping off the integrals with the parton distribu- 
tions functions. The remaining expression can be further simplified in the threshold region. 
For notational simplicity, we will perform these manipulations in the partonic cms frame 
where w = (1,0). We first shift the q integral by introducing the new integration variable 
q' = ( Xl P x + x 2 P 2 - 2Mw - q) = (y/§ - 2M)w - q and obtain 



(jppi 



^) = E E H *(ri J -^ Jl(Ew - q') Jdhe^W^(z,n). (2.62) 



i,i' 5=|s-s'| 



Here we defined the partonic centre-of-mass energy measured from threshold, E = Vs — 
2M = M(l — z) + 0(X 2 ). The function W contains soft fields that, by definition, vary 
significantly only on distances z ~ 1/A. Hence, only soft momenta (q' ~ A) contribute to 
the q' integration. On the other hand the function J, which is defined in terms of potential 
fields alone, is a function of q = (y§ — 2M)w — q' . In the partonic centre of mass frame we 
have by assumption X\P\ + a^-P? = Vsu> = pi + p 2 + k s = 2Mw + k' s with soft momenta 
k s , k' s , therefore q is also soft. Since potential fields with scalings (k ,k) ~ (A, vA) can 
depend on soft momenta only through their small time-like components, J is a function of 
go = Vs — 2M — q' alone. This enables us to perform the q integration in (j2.62p which 
sets z = in the argument of the soft function. We then obtain the final expression for 
the factorized short- distance cross sections: 



s+s' 





(Jpp/(s,/i) -- 


= E E*2 

S=|s-a'| i,i' 


wj 


duj2 J k(E- 


UJ 

~2 


we 


defined the Fourier transform of the soft function 








W*'(u,ii) = 


f dz 

1 4vr 


e^ z °/ 2 W^(z ,0 


u) 



W£?(u,n). (2.63) 



(2.64) 

The formula (j2.42p . with the partonic cross section (12.631) . establishes the factorization into 
collinear (the parton distribution functions), potential (the function J) and soft (the soft 
function W) contributions for heavy particles produced in an S'-wave, and constitutes the 
main theoretical result of this work. Some further comments on the structure and validity 
range of the factorization formula as well as a comparison to previous results are given in 
section 12.51 below. The simplification of the colour sum over i, %' due to the existence of a 
diagonal basis [37] is reviewed in section 13.11 
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2.5 Some comments on the factorization formula 

2.5.1 Gauge invar iance 

Let us briefly comment on the gauge invariance of the ingredients in the factorization 
formula. The hard function is defined in terms of on-shell fixed-order scattering amplitudes 
projected on a given colour and spin state. The gauge invariance of the hard function 
then follows from the gauge invariance of the on-shell scattering amplitude and the linear 
independence of the elements of the colour and spin bases. The effective-theory Lagrangian 
is invariant under the collinear gauge transformations of SCET [60J and the soft gauge 
transformations (I2.2T|) and f l2.20[) . Collinear gauge invariance is built into the formalism 
since the operators are constructed from the invariant fields C . The invariance of the 
components of the soft function (I2.64p under soft gauge transformations follows since the 
elements of the colour basis are invariant tensors satisfying 

Ji) _ Tj{R)\Tj{R')\ c (i) Tj(r) Tj(r') ( <y ggX 

C {a] — U a 3 b 3 U a A b A C {b} U b iai U b 2 a 2 - (Z.OOj 

Note that for this it is essential that only the soft function with vanishing spatial argument, 
W^, a (xo, 0, /z), enters the final factorization formula since for the soft function at arbitrary 
x the gauge transformations of the collinear, anticollinear and potential fields appear with 
different arguments x_, x + and Xq. Finally the collinear matrix elements and the potential 
function are defined in terms of the decoupled fields that transform trivially under soft 
gauge transformations. 

2.5.2 Subleading corrections to the cross section 

In our derivation of the factorization formula we relied on field redefinitions that decouple 
the soft gluons from the leading SCET and PNRQCD Lagrangians. The subleading effec- 
tive Lagrangians include higher-order potentials, but also interactions of the soft gluons to 
the potential fields via the x-E interaction mentioned above and analogous couplings to the 
collinear fields. These terms cannot be eliminated using the decoupling transformations. 
However, being sub-leading, these interactions can be treated as perturbations in (3 and, 
since the soft gluons decouple from the leading effective Lagrangian, the cross section can 
be written to all orders in the non-relativistic and SCET expansions in the schematic form 



a 



J2 H a [W a <g> J a ] , (2.66) 



with higher-order hard, soft and potential functions labeled by the index a. At NNLL 
accuracy, several effects could be relevant that are not included in the leading term a = 
in (I2.66P considered so far in this work. We now discuss these effects and show that they 
either do not contribute at NNLL or are incorporated in a straightforward way without 
spoiling the factorization. 

It is worth mentioning at this point that in standard situations of soft gluon resum- 
mation not requiring the consideration of Coulomb singularities, the expansion in powers 
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of 1 — z and a s can be considered separately. Thus the discussion of power-suppressed 
interactions is not necessary to any order in a s as long as one drops subleading terms in 
1 — z. In joint soft-gluon and Coulomb resummations the two expansions are linked by 
the counting a s / (3 ~ 1, so that 0((3) suppressed terms that would normally be referred to 
as power corrections now appear at the same NNLL order as a s terms, see (j2.7p . This is 
a complication characteristic of perturbative non-relativistic systems that have kinematic 
singularities which are stronger than logarithmic. 

Hard effects. Higher-order production operators such as the P-wave operator (12 . 18[) can 
appear at O(fl), which according to (12. 7p corresponds to a NNLL effect. These can include 
interactions with soft gluons through a spatial covariant derivative that are not removed 
by the decoupling transformation. However, since there is no interference between S- and 
P-wave production for the total cross section, these corrections are at least (D((3 2 ), beyond 
NNLL accuracy. 

Potential effects. Beyond the leading-order Coulomb potential in (12.191) . subleading effects 
lead to further potential interactions of the form 

d?r [^J b ] (f)5V abcd (r) WM] (0) • (2.67) 

The running of the strength of the Coulomb potential causes a NLL correction beginning 
with a 2 / (3 x log/3, which is accounted for by evaluating the coupling in the leading-order 
Coulomb potential at a scale of order m red /3, where m red denotes the reduced mass of the 
heavy-particle system, or by an explicit logarithmic correction to SV a bcd(r). Genuine loop 
corrections to the colour Coulomb potential lead to the substitution 

a.(l/r) _^ a^l/r) / + ^ + \ (2 6g) 



r r V 4ir 

in the Lagrangian (j2.19p . Only the one- loop correction Si is needed at NNLL, and con- 
tributes terms beginning with ot 2 s / '(3. In addition there exist the spin-dependent and inde- 
pendent non-Coulomb potentials of form a 2 /r 2 and (summarily) a s /r 3 , see e.g. [75], leading 
to NNLL terms beginning with a 2 log (3. Here the logarithm arises from non-relativistic 
factorization and is related to the fact that the short- distance functions H$, (/i) have poten- 
tial infrared divergences in addition to the familiar soft and collinear divergences. These 
additional potential terms in the PNRQCD Lagrangian do not involve soft gluon fields and 
the decoupling transformation can be applied as for the leading potential. Therefore these 
terms can be included in the evaluation of the potential function J, as was done in [15] to 
compute all log/3 and 1//3 enhanced terms of the it-production cross section at 0{a^). 

Subleading potential-soft interactions. We next discuss possible contributions to subleading 
terms in (I2.66P arising from the higher-order couplings of soft gluons to collinear and 
potential fields that are not decoupled by the field redefinition. In a diagrammatic language, 
these arise from the interference of subleading soft gluons coupling to the initial or final 
state with potential gluon exchange, that could contribute NNLL terms beginning with 
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a s /(3 x a s (3log ' (3. As shown in [45] using power-counting and rotational invariance, the 
fixed-order NNLO corrections of this form vanish for the total cross section. We will now 
use the effective theory language to show that the NNLL corrections arising from the 
subleading SCET and PNRQCD Lagrangians vanish to all orders in the strong coupling, 
so that the leading factorization formula (I2.63|) does not require modification at NNLL. 

We begin with the 0(j3) suppressed interactions in the subleading PNRQCD La- 
grangian [57l[70| [75]. 



(2.69) 



4nrqcd = ~9s [p\x)x ■ E s (x°,0)ij(x) + tfj'\x) x ■ E s (x°, 0)xf/(x) 
= -g s x-E a s {x\Q)J a {x) 
with the bilinear product of potential fields 

ja = ^\ T (R)a^ + ^/t T (K')<y . (2.70) 

We treat the chromoelectric vertex perturbatively, i.e. we consider contributions to the 
cross section where one of the matrix elements contains an insertion of the vertex x ■ E 
and evaluate these matrix elements with the LO PNRQCD Lagrangian. The first-order 
correction to (12.101) from the interaction (12.691) is given by 

AA^(pp'^HH'X)= C$ a} (») Jd'x(HH'X\T[C { X Q cB(^ M (mPP')EFT. 

(2-71) 
The expectation value is evaluated with the leading-order effective-theory Lagrangian, so 
the soft, collinear and potential fields decouple after the field redefinition. Under this 
redefinition, the chromoelectric interaction is transformed into 

ip j (x)x ■ E s (x ,0)ip(x) =ip m (x)x-£ s (x ,d)ip {0) (x) (2.72) 

with £ s = SljEgSw. Independent of the colour representation of the heavy particles, due to 
the identity (12.211) . which also holds for the Wilson lines, the components of the redefined 
chromoelectric fields are given by 

^(*o) = S® b E b s (x ) . (2.73) 

In analogy to ( 12.401) . the effective theory matrix element factorizes and the correction 
to the scattering amplitude can be written as 

AA^(pp' -► HH'X) = J2 cfafG") (O|0S 1Q » WtlJp') 

i 

d*xx(X\T[S:(x ,d)S§(0)}\0} (HH'\T[J^\x)^l^}\0) . (2.74) 

Repeating the steps leading to the factorization formula (12.621) and the subsequent discus- 
sion, we obtain an analogous formula with the replacement 
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^SwiSw) 1 J ^ {E - q0 ' k) I dz ° el9 ° z0 I d * x e ~ lk ' x s • ^ RaW {z °> x °^ ] + c - c • 

(2.75) 

Here we have defined the subleading soft function with an insertion of the chromoelectric 
field according to 

W2^ ( V,*°,A0 = P^Ty (0|T[4J ifc2} (^ )]Tg(x )4J fc3fc4} (0)]|0) (2.76) 

as well as a subleading potential function with an insertion of the bilinear f |2.70[) : 

J a { $ Ra (?. k) = P«? } j d'z e** f d'x e ik - x (0| [<) ai ^l 2 ](z)T[J<V {x) [^Vffl (0)] |0> . 

(2.77) 
Since the matrix elements are evaluated with the spin-independent leading order PNRQCD 
Lagrangian, the spin dependence simplifies in analogy to ( I2.47P as was used in (j2.75p . In 
the terms denoted by "c.c." in (I2.75P the operator S s is inserted in the anti-time-ordered 
product in the subleading soft function and the potential function is defined with the 
insertion of J^°> in an anti-time-ordered product with the annihilation fields. 

Because the subleading soft function does not depend on x, we can perform the integral 
over the spatial components of x in ( 12. 75ft and obtain the expression 



Iwr' si3)Ck) i r£>(E -"°' k ^ ' (278) 



which vanishes since (in the partonic cms frame) there is no three-momentum left that the 
integral over J^ can depend on. Hence we conclude that the correction to the cross section 
due to a single insertion of the NLO potential Lagrangian vanishes to all orders in the strong 
coupling constant. Similar to the diagrammatic argument in [35], an essential ingredient in 
the argument was the multipole expansion of soft fields when multiplying potential fields. 
Because of this, only E(x°) appears in the subleading PNRQCD interaction and, as a 
consequence, the subleading soft function depends only on x° which allowed to perform the 
simplification in the last step. The second ingredient is rotational invariance in combination 
with the independence of the potential function on any external three-momentum. 

Subleading collinear-soft interactions. The effects potentially relevant at NNLL arise from 
the effective Lagrangians at the next-to-leading order in the kinematic expansion. For the 
SCET Lagrangian of collinear quarks, the relevant interactions are given by [59JI60] 



4 1} = f {x^ W c gF; v W±) I i , (2.79) 

for n-collinear modes and similar terms involving transverse derivatives or factors of x± 
for the couplings to collinear gluons, and an interaction 

4« = Z^cWc q s + q s Wlilp ±c £ (2.80) 
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involving soft quarks but only collinear gluons. As argued already in [45], these vertices 
with transverse derivatives or x^_ cannot contribute to the total cross section since the 
initial-state momenta can be chosen to have zero transverse momentum. Then loop in- 
tegrals with transverse momentum factors vanish by rotational invariance in the plane 
transverse to the beam axis. Eq. f |2.80|) contains a vertex that describes the collinear split- 
ting of a quark into a gluon and a soft quark, which is also an 0(/3) term. However, a 
non-vanishing contribution to the squared amplitude requires two such vertices to contract 
the soft quark field, hence there is again no contribution at NNLL. Note that insertions 
of power-suppressed interaction Lagrangians lead to collinear matrix elements that do not 
take the form of the standard parton densities. At this point it appears that one must 
distinguish the collinear expansion in powers of k± ~ Ma/A ~ Mf3 from the standard 
collinear expansion in k± ~ Aqcd- However, as mentioned above, this issue is not relevant 
for NLL and NNLL resummations. 

To summarize, the factorization formula f 1 2 . 6 3 f) is valid at NNLL accuracy provided 
subleading Coulomb and non-Coulomb potentials are included in the computation of the 
potential function and the hard and soft functions are computed at NLO. 

2.5.3 Comparison to factorization at fixed invariant mass 

The result ( 12.631) can be compared to the analogous formula (11.11) given in [7J for the 
pair production of relativistic coloured particles with different four- velocity vectors W\,W2, 
where potential degrees of freedom do not appear. In this case the potential function is 
absent in the factorization formula and the soft function is given by (taking into account 
our convention for the antiparticles): 

S {ab} (z) = (OlT^^^.w^.wlWT^n,^^^^^ J(0)|0> . (2.81) 

In our case the extraction of the potential function J leads to a more complicated colour 
structure of the soft function W, but to a simpler kinematical structure due to the equal 
four-velocities of the heavy particles. Using the completeness relation of the projectors and 
the definition of the components of the soft function ( 12.561) . the limit of ( 12.811) for the case 
of equal four velocities, W\ = W2 = w, is formally related to the sum of the soft functions 
for the different final-state representations: 

S'tt / Li=i02 = 2_^ W u , a ■ (2.82) 

The resummation coefficients used in the Mellin space approach to threshold resummation 
have been calculated at the one-loop level for various colour representations [7J117J20] from 
the threshold limit of the soft function S. This amounts to taking the threshold limit after 
extracting the 1/e poles of the eikonal diagrams contributing to S. Applying this approach 
at the two-loop level [12] one encounters three-parton colour correlations proportional to 
log/3 that arise from an interference of soft and potential divergences, and lead to off- 
diagonal contributions in the colour basis constructed in [37] that diagonalizes the soft 
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functions W Ra . In contrast, the effective theory approach constructs an expansion in (3 
before taking the e — > limit, resulting in the tower of higher dimensional hard, soft and 
potential functions ( 12.661) . In this approach, the log /3-dependent divergences are attributed 
to higher- dimensional soft and potential functions due to the subleading x ■ E vertex in 
the PNRQCD Lagrangian and subleading potentials, while the leading soft function ( 12.56!) 
is diagonal to all orders. In order to resum all terms that are enhanced by powers of 
log/3 and inverse powers of (3 at a given accuracy one should use the formula (12.661) . 
while (II. ip is appropriate for partonic thresholds where the relative velocities of the final 
state particles are relativistic, as considered e.g. in the recent resummation of the invariant 
mass distribution of top-quark pairs in 



2.5.4 Soft-Coulomb factorization in Mellin moment space 

The enhancement due to single- Coulomb exchange has been included in the Mellin moment 
space approach to resummation as part of the hard function pj] , which implicitly assumes 
the factorization of soft and potential gluons. We now briefly discuss the relation of our 
factorization formula (12.631) to the Mellin space formalism and show that the convolution 
of the soft and potential function factorizes in Mellin space in the large N limit, which 
justifies the earlier treatment. We recall that the Mellin moments of the cross section in 
the variable z = 4M 2 /s are defined as: 

&(N,fi)= / dzz N - l a(s,fi). (2.83) 

Jo 

In order to Mellin-transform the factorized cross section (I2.63J) . we approximate E = 
VS — 2M « M(l — z) and make the transformation of variables li = 2M(1 —w). For stable 
particles and neglecting bound-state contributions, the potential function is non-vanishing 
only for E — u/2 « M(w — z) > so we can write 



&„,,($, fi)*2Mj2 H i> I ^ E J lW-^))W^(2M(l- W ),/;). (2-84) 



Here we have used z < w < 1 so that E — co/2 us M(w — z) ~ M{\ — z/w) and introduced 
a factor — rj 1. Since a convolution of the form f —f(w)g(z/w) factorizes under a Mellin 
transform, we conclude that, up to power suppressed terms, our factorization formula 
implies multiplicative soft-Coulomb factorization of the cross section in Mellin space: 

a pp/ (N, n) » Y, H %> E J k ( iV ) S u> (N, /i) , (2.85) 

S,i,i' R a 

where the soft function in Mellin space is given by 

/»1 /»oo 

S*r(N,(i) = 2M dzz N - 1 W* a {2M{l-z),ii)= dwe~^ u W^(u,fi) + O^- 1 ) , 
Jo Jo 

(2.86) 
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and 

JL(N) = [ dzz N - 1 J|JM(1 - z)) . (2.87) 

Jo 

This shows that previous treatments that put the one-loop Coulomb corrections into the 

hard function [11] can be extended to include the higher-order Coulomb corrections as well. 

2.5.5 Invariant mass distribution near threshold 

The combination of soft and Coulomb effects for the invariant mass distribution of top- 
quark and gluino pairs near threshold has been performed in several approximations as- 
suming soft-Coulomb factorization [23H2B]. We now show that a factorization formula for 
the differential cross section da / dMnw with M^ H , = (pi +P2) 2 , valid for M^h' near 2M, 
can be derived from our main result ( I2.63p . 
We first introduce the parton luminosity 

L pp >(r,fj)= dx 1 dx 2 S(xxx 2 -T) f p / Nl (x 1 ,/ji)f p '/ N2 (x2,ij) (2.88) 

Jo 

to express the hadron scattering cross section as 

a NlN2 ^HH'x = y^ / dTL pp >(T,n)a pp/ (Ts,n), (2.89) 

P ,p>= q , q , 9 J ^ 2 /° 

with a pp '(rs,fi) given by (12.63 p . Next we observe that the argument E — u/2 = y/rs — 
2M — to/2 of the Coulomb function J^ in that equation corresponds to the non-relativistic 
energy of the HH' pair in the partonic cms frame. In this frame the three momentum of 
the pair is soft and therefore makes a negligible contribution to the invariant mass, so the 
relation 

U — -o4N _ r— ,./o , zrVJ^fl^ 



M HH , = 2M + E-- + 0(M(3 A ) = V^ - u/2 + 0(M/3 4 ) (2.90) 

applies. We now change variables from 00 to Mhh 1 in (I2.63H and interchange the r and 
(implicit) 00 integration in ( I2.89P using 

dr (2.91) 

and the fact that for stable particles, and neglecting bound-state contributions, the Coulomb 
function has support only for positive values of its argument, which fixes the upper limit 
2E. This results in 

d<JN 1 N 2 ^HH'X V"^ f ATI \^pp'{ TS -.\A 
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dM HH > ,_ Jm 2 ,/s " ' dM H H' 

s+s' 

= E E E 2/ ^)E Jl(M HH ,-2M) 

p,p'=q,q,g S=\s—s'\ i,i' Ra 
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/ drL pp/ (T,fi)W l ^(2(V^-M HH/ ),fi) (2.92) 

JMl„,/ s 



with 
y '^ = £ S2^( A i)^Jl(M ra ,-2M)W«-(2(>/i-M OT ), A *), (2.93) 

HH ' S=\s-s'\ i,i> R a 

which is the desired result. This shows that Coulomb and soft effects factorize multi- 
plicatively in the invariant mass distribution near threshold, as assumed in [24T - f26] . Note 
that in the hadronic cross section only the soft function is averaged with the parton lu- 
minosity. Eqs. (12.921) . (12.93J) apply unmodified to bound-state contributions, in which 
case the invariant mass distribution extends below the nominal threshold 2M, and unsta- 
ble particles, where in addition the Coulomb function is evaluated at complex argument 
M HW -2M + i(T H + r H /)/2. 

3 Resummation of soft and Coulomb gluons 

The factorization formula (I2.63J) provides the basis for resummation of soft and Coulomb 
gluon effects in the soft function W and the potential function J, respectively. The soft 
gluon resummation is performed by solving evolution equations for the soft and hard func- 
tions [27H3T] . while the resummation of potential effects can be performed using results 
obtained for top-quark production in electron-positron collisions [761478] . In this section we 
provide the explicit results for the resummed soft and potential functions. In section 13.11 
we recall the colour basis constructed in [37] that diagonalizes the soft function to all or- 
ders in the strong coupling constant and quote the result for the one-loop soft function. 
In section 13.21 we relate the hard function to the colour- and spin decomposed Born cross 
section and the colour decomposed one-loop amplitudes at threshold to demonstrate that 
the implementation of the resummation formula requires only a standard calculation in 
fixed-order perturbation theory. In section 13731 we provide the leading order potential func- 
tion summing up multiple Coulomb gluon exchange, as required for the NLL resummation 
performed in section |H In section 13.41 we obtain evolution equations for the short- distance 
coefficients C and the soft function W that allow to resum soft-gluon effects. The anoma- 
lous dimensions required up to NNLL have already been collected in [37J. The final result 
for the resummed cross section using the momentum-space formalism of [3TH33] is pre- 
sented in section 13.51 with the explicit expression of the resummation exponent up to NLL 
relegated to appendix O Expansions of the resummed result to order a s , as required for 
the matching to a fixed-order calculation, are given in appendix [D] 

3.1 The soft function 

In the factorization formula (12.631) . the hard- and soft function are matrices in colour space 
in the basis of the tensors c® introduced in (12.141) . The soft function can be diagonalized 
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to all orders in the strong coupling for an arbitrary heavy coloured particle production 
process using the colour basis constructed in [57]. In this basis, the components of the soft 
function ( I2.56P can be expressed in terms of the soft function for a single heavy final state 
particle in the irreducible colour representation R a . This gives a precise meaning to the 
picture of soft gluon radiation resolving only the total colour charge of a heavy-particle 
pair at rest [11] . We review this construction here and quote the result of the one-loop soft 
function. 

It is useful to perform a decomposition of the product of the representations of the 
initial state and final state particles into irreducible representations: 



r®r' = ^r ai R® R' = ^R a . (3.1) 



It is intuitively clear that a final state pair in an irreducible colour representation R a can 
only be produced from an initial state system in an equivalent representation. In order to 
implement this picture formally, one forms pairs Pi = (r a , Rp) of equivalent initial and final 
state representations r a and Rp , treating multiple occurrences of equivalent representations 
in the decomposition as distinct, e.g. in the case of a symmetric or antisymmetric coupling 
of 8 (g> 8 — > 8. It has been shown in [37] that a basis respecting colour conservation ( 12.65)) 
can always been chosen by forming for every pair Pi of equivalent representations the 
associated basis element 

„(*) _ r«a ^R/3* (n r>\ 

y/djm{r a ) 

where the C^ are Clebsch-Gordan coefficients implementing a unitary basis transfor- 
mation from the basis vectors of the tensor product space R <S> R' to basis vectors of 
the irreducible representations R a , and analogously for the initial state representations. 
The Clebsch-Gordan coefficients, basis elements and projectors for all squark and gluino 
production processes at hadron colliders have been provided in appendix B of [37]. For 
squark-antisquark production from quark- ant iquark annihilation, the allowed pairs of rep- 
resentations are 

Pi G {(1,1), (8, 8)} (3.3) 

and the basis has been given already in (j2.16p . For the gluon fusion channel the allowed 
pairs of representations are 

PiE {(1,1), (8 5 , 8), (8 R , 8)} (3.4) 

and the basis is given in ( I4.13P below. 

Using properties of Wilson lines and the Clebsch-Gordan coefficients it was shown in [37] 
that the components of the soft function (12.561) in the basis ( 13.21) can be obtained from the 
soft function for the production of a single particle in the representation R a 

W& m {z,n) = (0|T[^ K 4 jfe2 ^ ij6 J(^)T[5 n , ail ^, a2J ^;l](0)|0) (3.5) 
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by contracting with appropriate Clebsch-Gordan coefficients 

«$"(<", A*) = /A . , \ A . , l 8 RaRfjl 8 Ra R, C[: a] W^ am {u^)C r ^ } . (3.6) 
y/ dim (r a ) dim [r a > ) 

As indicated, the elements are non-vanishing only if the irreducible representation R a is 
identical to both final state representations Rp and Rp/ in the pairs Pi = (r a , Rp) and 
Pi' = (r a ',Rpi) that define the tensors c^\ p- 1 '. This is intuitively plausible since we 
project on a specific final state representation so only production from initial states in an 
equivalent representation is possible. This shows that the soft function is automatically 
block diagonal in the basis (j3.2p where off-diagonal elements are only possible if several 
initial state representations r a are equivalent. For initial state quarks, antiquarks and 
gluons this only happens for two gluons in the initial state that can be combined into 
a symmetric and antisymmetric octet, see (j3.4p . Using Bose symmetry it is furthermore 
possible to show that the symmetric and antisymmetric octet production channels do not 
interfere so the soft function is diagonal, i.e. 

W*r(u,n) = Wf-(u,n) 8 iV 5 RaRr (3.7) 

The one-loop term in the loop expansion of the soft function, 

W*- ( W , n) = f; (^M) n Wt )R ° (*, /*) , (3.8) 

n=0 ^ ' 

was calculated in [37] and depends only on the quadratic Casimir operators C R of the 
representations of the initial state particles and the final state pair, in agreement with 
calculations for specific colour representations (TJUZUSO]. The result in position space can 
be written in terms of the variable L = 2 In ( lz °^ E ) and is given by 



wi 1)Ra {L,n) = (C r + C r r) ( - + -L + L 2 + ■— 1 + 2C Ra ( - + L + 2 1 , (3.9) 



where the normalization is such that W^ a (L,/i) = 1. For a final state singlet and a 
quark- ant iquark initial state this agrees with the Drell-Yan Wilson line [29], for a colour 
octet final state the result agrees with [21]. The Fourier transform of this result enters 
the factorization formula ( I2.63P and has been computed in [37]. Since in the momentum 
space formalism the solution of the renormalization group equation for the soft function 
is obtained using the Laplace transform of the momentum space soft function that can be 
obtained from the position space result by a simple replacement rule, the momentum space 
expression is not explicitly needed in this paper. 

3.2 The hard function 

In order to simplify the application of the factorization formula (I2.63|) we show how to 
bypass the matching to the effective theory used in our derivation and express the hard 
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function (12.601) directly in terms of hard-scattering amplitudes for the production of the 
heavy-particle pair. The potential and soft functions are universal functions, depending 
only on the colour quantum numbers. The only process dependent input for NLL resumma- 
tion are then the colour-separated leading-order production cross sections for the process 
of interest. For NNLL resummation, in addition the tree cross sections for individual spin 
states of the heavy-particle pair and the colour separated NLO amplitudes at threshold 
are required. We work in the colour basis that diagonalizes the soft function and denote 
the diagonal elements of the hard function by Hf . 

We recall that the S-wave matching coefficients CL'i are simply given by the scattering 
amplitudes at threshold, with polarization vectors removed and projected on a specific 
colour channel according to A{ a } = ^2 i CjKA^. Thus the scattering amplitude for a given 
colour configuration and fixed helicities and spins Aj can be written as 

,4(0 (pV* 9 -»- H^H'^) = 2 v /^^Cf4V) X^X&tXV . ( 3 - 10 ) 



where we have accounted for the normalization factors y/2mn in the definition of the non- 
relativistic states. Here the £ are the polarization spinors or vectors of the non-relativistic 
particles satisfying the completeness relation 



£)£#* = *<*■ (3-ii) 

A 

Computing the spin averaged partonic tree cross section for the production of a heavy- 
particle pair in a fixed colour state using the expression of the amplitude in the effective 
theory ( 13. lOf) we obtain 

°* [S) ~ 2s 8ns 4D r D r , ^ '^omlPP "> ^ )\ 

pol 

(3-12) 

In the second line we have approximated the prefactor using A(s, ra#, m 2 H ,) = (s—4M 2 )(s — 
{viih ~ m H') 2 ) ~ (s — 4M 2 )4m^m//' and the Cjji (//) should be computed in the tree 
approximation. Comparing to the definition of the hard function (I2.60J) . we obtain a 
simple relation between the leading order hard functions and the spin averaged total cross 
section for a given colour channel: 

The tree-level hard functions can therefore be obtained from the threshold limit of the 
Born cross sections in a specific colour channel. With the aim of extending the validity of 
the resummed expressions one may consider defining an improved hard function by using 
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full tree- level cross section instead of its threshold limit. In this way, some higher order 
terms in j3 are included, though not systematically. 

At NNLL spin-dependent hard functions are required in the tree approximation, which 
can be obtained by a formula analogous to (13. 13[) from the cross section for a fixed colour 
and spin channel. In the framework of a standard computation of the leading order cross 
section, the projection on the final-state spin can be performed introducing scattering 
amplitudes for a fixed final state spin S that can be obtained from the helicity amplitudes 
for the production of the HH' pair according to 

A (i \p Xl p ,X2 -> (HH') S ' X ) = J2 A^AIAgA^^W* 2 -► H Xi H' X4 ). (3.14) 

A3A4 

The iVs^AIAsA^ are Clebsch-Gordan coefficients that can be found for the case of spin 
one-half particles, e.g. in [75] . 

Furthermore, at NNLL accuracy also the one-loop hard function is required. Since the 
spin-dependence of the potential function is already an 0(a 2 s ) effect, one may use the spin 
summed one-loop hard functions H\ at NNLL without formal loss of accuracy. The one- 
loop hard function is given by the interference of the Born and the one-loop amplitudes 
for a given colour channel, evaluated directly at threshold: 



#fV) = ^\K 1 ,A 2 ^ic^^^c^*M) 



2s L 



SmuTriH'S AD., _ , 

pol 



^- £ 2 Re (4t>P -► HH')A^ ( PP ' -► HH')) . (3.15) 



Here -4nlo is the UV-renormalized one-loop amplitude evaluated directly at threshold 
with IR singularities regularized dimensionally and subtracted in the MS-scheme. This 
is therefore the only process independent input required for NNLL resummations and is 
far simpler to compute than the full NLO cross section. Alternatively, the hard function 
can be extracted from the threshold expansion of the NLO partonic cross sections in each 
colour channel [15] . These are available for ti production at hadron colliders [80] but not 
yet for the production of squarks or gluinos. 

3.3 The potential function 

We now discuss the relation of the potential function J to the Coulomb Green function 
familiar from PNRQCD. The momentum-space potential function (12.591) can be written in 
terms of the correlation function 

pol, colour 

*ze** (0|[<l^l](z)[^3< o r 2 ](0)|0) , (3.16) 
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where the matrix elements are evaluated with the leading order PNRQCD Lagrangianjj 
In the second line we used that in this approximation the soft gluon fields are decoupled 
from the fields ip^ so we can replace the sum over the two-particle Hilbert space by a 
sum over the full Hilbert space and use the completeness relation. The definition of J 
given here is sufficient for 5-wave production of the heavy particles; for P-wave production 
also expectation values of fields with derivatives have to be considered. As explained in 
section 12.5.2] effects of the x ■ E vertex in the sub-leading PNRQCD Lagrangian are not 
included in the calculation of the potential function J at higher orders but lead to additional 
contributions to the factorization formula with generalized soft and potential functions, to 
be calculated with the leading effective Lagrangians. As shown in section [27572] these terms 
are only relevant for the total cross section beyond the NNLL order. 
Defining a tensor product notation for the decoupled potential fields, 

&®iP') aia2 (t,r,R) = ^(t,R + 0V4 O) M~ f) , ( 3 - 17 ) 

with the relative coordinates r and the cms coordinates R and using the projectors onto 
the irreducible representations R a introduced in (12.461) . we can perform the colour decom- 
position 

Ra Ra 



This is analogous to the singlet-octet decomposition discussed in [67] (see also [81] . in 
particular eq. (48)). With this notation, the colour structure of the Coulomb potential 
simplifies to 



Ha, Ftp 
Ra 

(3-19) 
Here the coefficients of the Coulomb potential are defined by the relation 

TXT£fP^ 2a3a4 = D Ra P^ (3.20) 

and we have used the projection property 

pRa pR/3 __ X pRa (O n-l\ 

a\a,2b 1 b2 b\b2C\C2 u R a Rfi r \\a2c\c2 " yo.z,L) 

Using the simplification of the Coulomb potential, the leading order PNRQCD La- 
grangian for the decoupled fields can be written in the tensor-product notation as a sum 



9 As discussed in section 12.51 to reach NNLL accuracy, the LO Lagrangian is to be supplemented by 
NLO Coulomb and the leading non-Coulomb potentials. 
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over the irreducible colour representations, 

Ha V red K ) / ^ 22 j 

+ fd 3 r[{iP®iP') RaJ! (t,-r/'2,0)] (^^-) [(?p®?p') Ra (t,f/2,0)] 

with the reduced mass m red = m H mH' / (% + %')• The kinetic term can be written in 
terms of the tensor field as shown here after a projection on the two-particle sector of the 
theory [67j[8T]. 

For top-antitop and squark-antisquark production the irreducible representations are 
given by3(g)3 = l©8 and the coefficients are familiar from quarkonium physics (recall 
our sign convention for the Coulomb potential): 

N 2 - 1 
D\ = —Cp = ^— — (attractive), 

2N C 



D, 



^F — 



2N< 



(repulsive) . (3.23) 



c 



The explicit values of the coefficients for all remaining representations relevant for the 
production of coloured SUSY particles are collected in appendix iBl 

Since the Lagrangian is now diagonal in colour and spin, the leading-order Coulomb 
Green function is of the form 



K'iJ - 

R, 

with the correlation function 



J W(?) = E P W *«!«.*»«. jRa (^) ( 3 - 24 ) 



J R «(q)= /rf 4 ze^(0|[^ ,(0 V (0) ](^)[^ (0) V (0)t ](0)|0) . (3.25) 



In this expression the fields carry no colour and spin indices any more and the correlation 
functions are to be calculated with the term in fj3.22j) corresponding to the representation 
R a with Coulomb potential a s Dn a /r. Since the annihilation fields ip annihilate the vacuum, 
we can replace (^//)(z)(^t^ ^)(0) = [(ipip')(z), (ip^ip t)(0)] in the vacuum expectation value 
and express the correlation function in terms of the imaginary part of a time ordered 
product 

J R «( q ) = /d 4 ^e^2Im(0|T[V'' (0) V' (0) ](^[V' (0)t V'' (0)t ](0)|0) 

J (3.26) 

= 2ImG Ra(0 \o,0;E). 

Here we introduced the zero-distance Coulomb Green function of the Schrodinger operator 
— V 2 /(2m re d) — (—D Ra )a s /r, i.e. the Green function G c a '(r^r^E) evaluated at fi = 
r 2 = 0. 
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Using the representation of the Green function given in [82], the MS-subtracted zero- 
distance Coulomb Green function including all-order gluon exchange, reads as follows |7U] : 



yR a (0) 



G£ (0,0;£) 




(3.27) 



Here 7^ is the Euler-Mascheroni constant and one should apply the prescription E — > E+iS 
for stable heavy particles and E — >• E + i(T H + T H /)/2 for unstable ones. 

For positive values of E the potential function evaluates to the Sommerfeld factor 



(2m red ) 2 7iD Ra a s 



-1 



J Ra (E) = ± ^_ f!«nf e ^«-«.V-^-lj E>0. (3.28) 

If the potential is attractive, D Ra < 0, there is a sum of bound states below threshold 
given by 

J h £ md (E) = 2JT$(E- E n ) f 2m rodhD Ra )a s \ E<Q ^^ 

n=l \ n J 

with bound-state energies 

2m red a 2 D 2 R 
E n = ! ^ 5a - (3-30) 

A series representation of the imaginary part of the Coulomb Green function for finite 
widths and arbitrary energies can be found in |23j . 

The above results suffice for resummation with NLL accuracy as performed in section HI 
For NNLL accuracy in the counting (12. 7p one has to resum in addition all a s x (ct s //3) n 
corrections as well as the non-relativistic logarithms of the form a 2 \n(3, a^ln (3, . . . . The 
fixed-order NNLO corrections of this form have been obtained in |45j. An analytical 
result resumming all a s x (a s /(3) n terms was obtained in [76] and is quoted e.g. in [25] . 
Resummation of the non-relativistic logarithms requires the generalization of results for 
top-quark pair production obtained e.g. in [83j[8l] to arbitrary colour representations. 

3.4 Evolution equations of hard and soft functions 

In the momentum-space approach to threshold resummation one calculates the short- 
distance coefficients C®(/jl) at a hard scale fih ~ 2M and the soft function W(u,f/,) at 
a soft scale of the order of /j, s ~ MP 2 . Threshold logarithms log/3 are resummed by using 
evolution equations to evolve both functions to an intermediate factorization scale fif. In 
this subsection we will provide these evolution equations in the colour basis ( 13. 2 p that diag- 
onalizes the soft function. In this case the structure of the evolution equations is similar to 
those in the Drell-Yan process [33, ,35] and the complications of matrix-valued anomalous 
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dimensions depending on the kinematics [TIGS] are avoided. The resummed cross section 
obtained from solving the equations is given in 13.51 

The evolution equation of the hard coefficients has been obtained in [37] from the results 
of [38] for the IR structure of general massive QCD amplitudes: 



dliafi 



{<*} 



(X'r + C/) ( hi ( — r ) ■-,!<) +/'7rC' / ,.„ 



V \ Mi) 



+7nqi>)- 



'" ~ 1 2 ' cusp ^ rTUr 'M m l — ~ ) ~ "" ) ^ '" w «o t /. l ^j^lf 

(3.31) 

Here Casimir scaling was used to write the cusp anomalous dimension for a massless parton 
in the representation r in the form r£ usp = C r 7 cusp which holds at least up to the three-loop 
order. At least up to the two-loop level the anomalous dimension jY can be written in 
terms of single-particle anomalous dimensions: 

it = i r + i J + it- ( 3 - 32 ) 

The anomalous dimension ^ H a s is related to a massive particle in the final state representa- 
tion R a in the pair Pi = (r' a , R a ) defining the colour basis element & 1 ' with index i. While 
the one- and two-loop anomalous-dimension coefficients Y of massless quarks and gluons 
have been available for some time, the two-loop results for the heavy particle soft anoma- 
lous dimension j H a s have become available only recently [3TPXJ . The one- loop coefficients 

of the cusp and soft anomalous dimensions are simply 7cus P = 4 and 7^™ = —2Cn a . 
The anomalous dimensions Y related to the light partons do not obey Casimir scaling 
already at one loop. The explicit one- and two-loop results for all remaining anomalous 
dimensions appearing in this section are available in the literature and have been collected 
in [37]. We observe that, as noted for the production of a colour octet final state particle 
in [21] , the imaginary part in the evolution equation (13.311) can not simply absorbed in the 
argument of the logarithm by the continuation M 2 — > — (M 2 + iO) which complicates the 
resummation of "7r 2 -enhanced" terms compared to colour singlet final states as in Higgs 
production [3UES]. We will not consider such a resummation of constant terms here. The 
evolution equation of the hard-functions is obtained from (I3.3ip as 

dJ^/O*) = (WCV + C r> ) In (^Pj + 2 7 f) HfQi). (3.33) 

As discussed in [37] and section I2.5.2[ starting from the two-loop order there are IR 
divergent contributions to the short- distance coefficients that are related to UV divergences 
of insertions of the non-Coulomb PNRQCD potentials in the extended factorization for- 
mula f)2.66p . These divergences cause additional scale dependence, not included in (J3.33J) . 
which is cancelled by a non-trivial scale-dependence of the potential function including the 
non-Coulomb potential insertions. The factorization scale dependence in the separation 
of H and J is related to additional non-relativistic log /3 terms of the NNLL order. Since 
in this paper we do not consider the resummation of these non-relativistic logarithms, in 
order to obtain the evolution equation of the soft function, these contributions to the scale 
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dependence of the hard function have to be dropped as done in (13.331) . Beyond NNLL fur- 
ther complications can arise from the structure of the extended factorization formula ( I2.66P 
including terms with higher- dimensional soft functions, but the discussion of resummation 
beyond NNLL is beyond the scope of this paper. 

The evolution equation of the soft function can be obtained from that of the hard 
function using the factorization scale independence of the hadronic cross section. Consis- 
tent with our discussion above, we consider the potential function to be scale independent 
which is appropriate for the NLO potential function quoted in [25J. Scale invariance of the 
total cross section and the known factorization scale dependence of the PDFs implies the 
evolution equation of the partonic cross section 

d s~^ f l dx 

' z ^) = ~z2 —(P p /p(x) + P p// p / (x))a P p / (z/x,fi), (3.34) 



-O rm > 



dha.fi ' 4rf J z x 



p,p 



where P p /p (x) are the Altarelli-Parisi splitting functions for the splitting of a parton p into 
a parton p. Consistent with the derivation of the factorization formula in section [2] from 
the on-shell scattering process pp' — > HH'X at production threshold, we use the x — > 1 
limit of the splitting functions for a parton p in the colour representation r, 

P p/p (x) = (2T r cmp (a s )j^r- + 2j^(a s )S(l - x)\ 5 PP + 0(1 - x) . (3.35) 

As discussed in section I2.5.2[ subleading collinear terms could potentially be enhanced 
by the Coulomb singularity, so care in dropping higher order corrections in equations like 
this is required. In the present case, however, the corrections to f !3.35[) are of the order 
1 — z = (3 2 so they they are not relevant at NNLL accuracy. We recall the property of the 
plus distribution, 



dx f(x) 



x 



1 JV~\ ffi\ l-z 



/(*)-/a) r^ /(i) 



dx ^T^x) — h dx T^x)- (3 ' 36) 



Inserting the factorized form of the partonic cross section into ( 13. 34ft and making use 
of the evolution equation of the short- distance coefficients ( 13.31D leads to an evolution 
equation of the soft function. Using the relations s ~ AM 2 , (1 — z/x) ~ (1 — z) — (1 — x) 
valid in the x — > 1, z — > 1 limits and introducing the variable u = 4M(1 — x), the 
resulting equation can be written in the same form as the equation found for the Drell-Yan 
process 



'' W*"(u;,fi) = -2 



J - .'us,, " ^„)l"l^)+27fe 



W«(u,fi) 



-2(r_ + r Lp )/ ^ w^^-wf-M (3 _ 37) 



din ft 

- 2 (rL, P + r^)/"^5£^_ w , 

with the anomalous- dimension coefficient 

7fe = 7f + 7 ' r + 7^'- (3.38) 
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These results hold for the colour basis (13.21) that diagonalizes the soft function to all 
orders [57], where only the diagonal elements Wj = W„ of the soft function have to be 
considered. Analogously to the anomalous dimension jf of the hard function (13.32 p . at 
least up to the two-loop level the anomalous- dimension coefficient of the soft function f l3.38|) 
can be written in terms of separate single-particle contributions 

lw7 = it + Tl + Is (3-39) 

with 

7 I = 7 r + 7* r - (3.40) 

The anomalous- dimension coefficients 7^ vanish at one-loop level so that 7^" = 7^™ = 
—2Cn a . Taking the Fourier transform of the evolution equation of the position space soft 
function given in (3.30) of [37] one obtains the momentum-space evolution equation quoted 
in ( 13.371) . In order to obtain this result, the logarithmic term multiplied by r cusp has to be 
treated with care, for instance by writing lnx = lim^o -(1 — x~ e ) (see also [28J). 

The terms proportional to the cusp anomalous dimensions in the evolution equations 
are related to the resummation of double logarithms In (3, while the coefficients jf and 
7 W ™ are related to single logarithms. For LL summation one needs the one-loop cusp 
anomalous dimension while all other quantities enter at LO. For an NLL resummation 
the required ingredients are the two-loop cusp anomalous dimensions and the one-loop 
anomalous dimensions 7^ and 7 W ". For a NNLL resummation one needs in addition to 
the three-loop cusp anomalous dimensions and the two-loop anomalous dimensions, the 
one-loop soft function (13. 9p and the one- loop hard function (I3.15p . 

3.5 Resummed cross section 



As mentioned in the beginning of section I3.4[ the resummation of soft-gluon corrections in 
the approach of [3"TH3"5] is performed by calculating the hard and soft functions at scales \ih 
and \i s that minimize the radiative corrections to these quantities and subsequently using 
the renormalization group equations (13.331) and (I3.37P to evolve to a common scale \if 
to compute the partonic cross section (I2.63P and perform the convolution with the PDFs 
evaluated at the same scale (see figured]). 

Since we have determined a colour basis that diagonalizes the soft function to all orders 
in perturbation theory for all hadron collider processes of interest, the evolution equations 
have the same form as those for the Drell-Yan process [33]. The resummed hard function 
solving the evolution equation f!3.33|) is given by 

(4M 2 \ ~ 2ar ^ h '^ 
— H?(jjih) , (3.41) 

with the functions S, af and ar defined as [33] 

~ a * M , ^ usp (a s ) + r msp (a s ) f"° da' s 



S(/j, h ,jj,) = - da 

J a,(u.h) 



own,) 2/3(« s ) J*.Q* h )P(<x a 
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W^fat*. 
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Vs 



Figure 1: Sketch of the resummation of soft gluon corrections using RGEs. 



«r(/ih,/i) 



as ^ T r (a ) + T r ' (a 



(3.42) 



Here a s (/i) represents the QCD coupling in the MS scheme and (3(a s ) the corresponding 
(3- function. Explicit results for the functions (3, ar and the Sudakov exponent S up to the 
NLL order as needed for section @] are collected in appendix E2 expressions up to the N 3 LL 
order can be found in [33] . 

The evolution equation of the soft function (13.371) can be solved in momentum space 
by introducing the Laplace-transform with respect to the variable s = l/(e 7B /ie p / 2 ), 



°i 



(P,V) 



due 



-W/f>,/i) 



(3.43) 



F^Q 



where we have defined the MS-renormalized soft function W{" . In practice, it is not 
necessary to perform the Laplace transform explicitly, since the soft function in position 
space W i a (z ,p) depends on the arguments solely through the combination 



izou.e'fz _ L/2 



(3.44) 



It is then easy to see that the function s { a (p) is obtained by simply replacing L — > —p in 
the MS-renormalized result for the soft function in position space. From the results of the 
tree-level and one-loop soft function in the diagonal basis (13.21) quoted in (I3.9P we obtain, 
discarding the 1/e poles due to the MS subtraction, 



~{0)R a 



(/>,/*) = 1. 



S^ (P, V) = (C r + C r> ) (p 2 + y)~ 2Cr « (P - 2 ) • 



(3.45) 
(3.46) 
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The resummed soft function can be written in terms of the Laplace transform with the 
argument p replaced by the derivative with respect to an auxiliary variable rj that is set to 
r] = 2ar(/i s ,/i) after performing the derivatives: 



Ra 



\zRc 



W^ es (u, fj.) = exp[-4S(/z fl , /i) + 2a^(> s , [i))sf<* (d v , [i 



1 



2 n 



9{u)- 



-27 B 77 



(3.47) 



Here the functions a v ^ i are defined in analogy to the functions (I3.42p with the obvious 
replacement of the anomalous dimensions jf — >■ 7^,". For rj < the solution (I3.47P should 
be understood in the distributional sense [32J as will be discussed below. The evolution 
equation of the soft function (13.371) can also be solved in Mellin moment space [29] and the 
standard resummation formula in Mellin space is reproduced. The precise relation between 
the quantities in the momentum space approach to resummation and the Mellin moment 
space formalism is discussed in [3"2" | l3"3" t l3"5 ] l3"7]. 

Using the solutions of the RGEs to evolve the hard and soft function in the factorized 
expression for the partonic cross section (I2.63P to the common factorization scale [if we 
obtain the result for the resummed cross section 



■ res 1 
PP' 



(T- f {S,[lf) 



Yl Yl H ?^h) u^m, [i h , [t s , [if) 



S=\s-s'\ 



2M 
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-2 n 



xsf-{d m n s 



-ZlEV 



T(2 V ) 



duj 



JL(e- 







U) 



[Is 



2,, 



s+s' 



J2 ^iJfw^fM,^^,^ 

5=|s-s'| i 

f°° Ji (E - %) 
Jo w 



-) *?■ 



2M 



2 In 



d v ,[i t 



e -2-y E rj 



(3.48) 



where have defined the evolution function 



Ui(M,[i h ,/j,f,/j, s 



f'l 



3)'"- 



A{S([i h ,[if) - S([x s ,[x f )) 



2af(/i h ,/i s ) + 2a' p ' r ([x s ,[Xf) + 2a </>,r '(/i s ,/i / ) 



(3.49) 



and now rj = 2ar(/i s ,/i/). The sum over the final state representations R a in the factor- 
ization formula (12.63}) has disappeared in the colour basis (13. 2p since there is a unique 
final state representation for each term in the sum over i (see the expression (13. 7p for the 
diagonal soft function). The evolution function Ui has been simplified using the identity 



a T ([th, Hf) + ar(yU/, Ps) = ar(Hh, Us) 
It could be further simplified using the identity 



S(Hh, Vf) - S{fi s ,fif) = S([i h , [i s ] - a r ([i s , [if) In 



til 

fls 



(3.50) 



(3.51) 
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to obtain 



Ui(M,fi h ,fif,fi s 



AM 2 



2\ -2a r (Mh,/i s ) 



x exp 



AS(nh,p. s ) 



2aY(ii h , n s ) + 2a^(/i s , N ) + 2a*' r# (// s , /*/)! , (3.52) 



which makes it clear that the \x$ dependence is related to the parton density functions. 
Eq. (13.511) is formally satisfied at whatever N fc LL accuracy employed, but not strictly valid 
in an expansion in the strong coupling since S and ar are evaluated at different orders. In 
the numerical evaluation we do not use this simplification but use (13. 49ft . 

A similar convolution of the resummed Coulomb Green function as in (I3.48P was en- 
countered in [61J where the mixed Coulomb and one-loop soft corrections to W-pair pro- 
duction were considered. The formula (I3.48J) generalizes the result in [61] to QCD and 
by resumming soft-gluons to all orders. An analogous formula was already used in [M] to 
resum QED corrections of the form Ya(Yw/Mw) for iy-pair production in electron-positron 
collisions. 

At the NLL order only the trivial, LO 0(a®) soft function, the spin-independent LO 
potential function and the leading-order hard functions are required so that the resummed 
cross section is given by 






Y,H\ Q) {n h )Ui(M,n h ,n a ,n f ) 



-2lEV 



r(2ij) 



00 J R (E- %) 
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-^)\ (3.53) 
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As mentioned before, for 77 < the kernel - (^) should be understood in the distri- 
bution sense [32]. Neglecting heavy particle decay and bound-state effects, the potential 
function J^ a vanishes for E < and the u integral is defined as the star-distribution 
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(3.54) 



where the double subtraction is sufficient to render the integral convergent for 77 > — 1. 
The poles of the original integrand at 77 = and t\ — —\ show up explicitly in the second 
line but are cancelled by the overall prefactor 1/T(2r]) in (I3.48P and (I3.53p . 

Note that the all-order solution ( I3.48P to the RGEs does not depend on the hard and soft 
scales fit and fi s , but truncating the perturbative expansion at a finite order in a s introduces 
a residual dependence on these scales, which is of higher order in a s . As pointed out 
the standard resummation formula in Mellin space following [l]|2] corresponds to 



m 



the implicit scale choices \ih = \if and \i s = M/N for the iV-th Mellin moment. The explicit 
appearance of these scales in the momentum-space formalism allows to use the residual 
scale dependence of the resummed cross section to estimate the remaining uncertainties 
from uncalculated higher order corrections. Our choices for these scales and the estimates 
for the remaining scale uncertainties are discussed in section 14.21 
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3.6 Matching to the fixed-order NLO calculation 

For the implementation of the NLL resummed cross section, we match to the fixed-order 
NLO calculation valid outside the threshold region. The matched result for the partonic 
cross section is then given by 



~ matched / %\ 



^(*)-C L(1) (s) 



6(A-[V~s-2M}) + a^°(s 



Aa^(s,A) + a^°(s) 



(3.55) 



where a??}? (s) is the fixed-order NLO cross section obtained in standard perturbation 



theory and a , is the resummed cross section expanded to NLO. The result is given 
explicitly in ( ID. II) in appendix [D] Since the resummed NLL cross section is expected to 
be a good approximation to the total cross section only near the threshold, we allowed for 
a cutoff E = Vs — 2M < A to switch off the resummation outside the threshold region, as 
was done in |14j . In the calculation of squark-antisquark production we choose, however, 
not to introduce this cut-off as discussed in section 14.21 

The total hadronic cross section at NLL is then obtained by convoluting f)3.55|) with 
the parton luminosity, 



<ttfmW= E f dTL pp ,(T, H )a™ tched (sT, H ), (3.56) 

P , P >= q , q , 9 J ^ 2 I° 

where the parton luminosity is defined in terms of the PDFs in (12. 88 p . 



4 NLL resummation for squark-antisquark produc- 
tion 

In this section we use the results of section |3]to perform the NLL resummation for squark- 
antisquark production in proton-proton and proton-antiproton collisions. According to 
the systematics introduced in (j2.7p . NLL accuracy requires the combined resummation of 
Coulomb gluons and soft gluons. Therefore our results extend previous ones on squark- 
antisquark production at higher order [171420] . that treated soft and Coulomb resummations 
separately or used a fixed-order expansion. In section 14.11 we perform the matching to the 
effective theory, and compare the threshold approximation of the LO and NLO cross sec- 
tions to the full results. In 14. 21 we discuss our choices for the hard, soft and Coulomb scales 
appearing in the resummed cross section and present predictions for squark-antisquark 
production at the Tevatron and the LHC at cms energies of 7, 10 and 14 TeV. We also 
compare to results from soft gluon resummation in the Mellin-moment approach [20] . 

4.1 Threshold expansion of the squark-antisquark cross section 

In the following we perform the matching from the MSSM to an effective theory containing 
only (anti-)collinear partons and non-relativistic squarks and antisquarks, and derive the 
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hard functions for the partonic subprocesses. As discussed in section [3T2l the hard functions 
required for NLL resummation can be inferred from the leading-order cross sections for the 
colour- singlet and octet production channels [TH|n^]. Therefore the explicit matching to 
the EFT is not strictly necessary, but is included here in some detail in order to provide an 
illustration for the somewhat abstract discussion in section |2j We also use the general result 
for the threshold expansion at NLO (ID.3J) to reproduce the known threshold behaviour of 



the NLO cross section [55] . 

We consider squark-antisquark production at hadron colliders, which at leading order 
proceeds through quark-antiquark and gluon-induced subprocesses: 



g(h)g(h) -> qaApi)Qa 2 j(P2) 



(4.i; 



where i,j,k,l denote the quark/squark flavours and aip = L/R label the scalar partners 
of the left /right-handed quarks. The relevant Feynman diagrams are shown in figure [2j 
We follow the setup of [85] and take all squarks to be mass degenerate and do not consider 
top squarks. The first restriction is done for simplicity (and was also used in other recent 
works on higher-order soft-gluon effects [TTH2U] ) and is not essential in our formalism, since 
the results in section [3] include the case of different final state masses. The application to 
top squarks would require an extension of our framework, since in the quark-antiquark 
production channel they are produced in a P-wave [86J. Resummation for top squarks was 
performed recently in [87j assuming that the resummation formalism for S- wave production 
can be applied to P-wave dominated processes. Since, in the present work, we are concerned 
with the total cross section, there are no sizable corrections from the finite squark width 
T as long as T <C M unlike the case of the invariant mass distribution near threshold. 
Contributions to the cross section from below the nominal production threshold will be 
included through the bound-state contribution to the Coulomb Green function. From a 
recent study of threshold effects in gluino production J2S] we expect additional finite- width 
effects at most of the order of these bound state corrections. 

To obtain the leading-order production operators we evaluate the tree-level scattering 
amplitudes for the processes ( 14. ip at the production threshold, (k% + k 2 ) 2 = ^m 2 ^ and 
determine the short- distance coefficients from the matching condition (I2.17P in order to 
compute the hard functions ( I2.60p . 

4.1.1 Production from quark-antiquark fusion 

The tree-level scattering amplitude for quark-antiquark induced squark-antisquark pro- 
duction at threshold is dominated by t-channel gluino exchange, if quark and squark (and 
antiquark and antisquark) flavours are identical. Near threshold the full-theory tree- level 
scattering amplitude for the quark-antiquark channel is 

^4{°a}(<Mj -> qkLQm) = - —TT -^ S ikSjiT^ a4 T^ ai v{mgn)P L u{rrign) . (4.2) 

m q i" m g 
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Figure 2: Leading-order Feynman diagrams for the quark- ant iquark (top) and gluon-fusion 
(bottom) induced production of squark-antisquark pairs. 



Here we use the usual projectors Pr/l = (1 ±7 5 )/2. There are analogous expressions with 
left and right labels exchanged. 

In the effective theory we introduce fields ipi (^) that annihilate a left squark (right 
antisquark) of flavour i. The quark fields are given by the SCET fields (j2.25p . We can 
then match the amplitude onto the effective theory according to the condition ( j2.17p . with 
the production operator given by (no sum over flavours is implied) 



e> (0) 



l(hcwWc) a2 (w&, c;ai ) ai ] U\, a ^% 



j;a 4 



(4.3) 



Here £j jC and ^ are the SCET fields describing collinear and anticollinear quarks with 
flavour i. From the matching condition (12.171) we can read off the matching coefficient, 
taking the non-relativistic normalization factor 2iriq into account: 



a 



qq,{a;a} 



9s m g 



m q{. m \ + m l) 



0,20,4 0301 K^Ljotzai 



(4.4) 



We now introduce the decomposition of the coefficients into a colour basis (j2.14|) . but 
leave the spin indices open. For the case at hand, the colour basis has been given already 
in f)2.16p . We find for the matching coefficients for the two colour channels: 



C 



(i) 



-c f 



Ana s m g 



rriq(m~ + m|) 



L)aia\ 



c 



(2) 



C F 47ia s m~ g 
2N C mdml + m? N 



> 2^ V" 1 L)a 2 a 1 

r mqyiibg ■"■'■ 



(4.5) 



From the definition (I2.60P and the polarization sum flA.Tj) we obtain the diagonal elements 
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of the hard function 10 !. 
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\m\ + m?) 



(4.6) 



(4.7) 



for the two colour channels. From ( 13.131) we obtain the total partonic Born cross section 
at threshold 

M I _ ^fP 

..- -~ \s=4m 
><liLqjR 



a 



iiij- 



27iN r 



na s m g 
[m\ + mf, 



(4- 



Adding the same expression with left and right squark labels exchanged reproduces the 
result given e.g. in [85]. In figure |3] we compare the partonic cross section at threshold to 
the full partonic LO cross section [85], for a squark mass of 1 TeV at the LHC with 14 TeV 
cms energy. We plot the integrand of the convolution in the formula for the hadronic cross 
section (j2.89p . i.e. the product of the partonic cross section and the parton luminosity, as 
a function of f3, taking the Jacobian |^ into account: 



da 



(0) 

d/3 



iPm 2 . 



s(l-/3< 



—L qq (P,n f ] 



(o) 



a 



(4.9) 



In the plot we use the MSTW2008LO PDFs [SS] with /x / = m~ q . It can be seen that the 
threshold approximation is applicable for /3 < 0.3, while the main contribution to the cross 
section comes from the region /3 rj 0.6. The results shown here are for a gluino-squark 
mass ratio of 1.25. For other gluino masses the discrepancies in the region /3 > 0.3 can be 
somewhat bigger. 

An approximation to the NLO cross section can be obtained by expanding the re- 
summed cross section ( I3.48P to order a^. The resulting expression can be written in the 



r (i) 



Mas f(l) 



form Acr , = crL/f 1 /™/ + O(P), where the scaling functions frj for an arbitrary colour 



(i) 



w 



)>)>' 



channel are given explicitly in (ID.3[) . Summing up the two colour channels we obtain the 
threshold expansion of the scaling function for the quark- ant iquark channel 
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(4.10) 



The logarithmically enhanced terms and the Coulomb correction agree with the results 
of [85]. The one-loop hard coefficient h^ is currently unknown in analytical form. In the 
right-hand plot in figure [3] we plot the full NLO corrections to the partonic cross section 
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For simplicity, in (|4.6[) and f|4.15[) we set l/(2s) ~ l/(8m|). However, in our numerical implementation 
of the resummed cross section (|3.53[) this pref actor is kept unexpanded. 
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Figure 3: Partonic cross sections for q^ — > %%, multiplied with the quark-antiquark 
luminosity and summed over quark and squark flavours, for rriq = 1 TeV and rrig = 
1.25 TeV at the 14 TeV LHC Left panel: exact leading order result (black/solid) and 
threshold approximation (blue/dashed). Right panel: exact NLO result (black/solid) and 
the approximation based on f)4.10p with h q <j{lJ) = 0, using the threshold approximation of 
the tree (red/dot-dashed) and the full tree (blue/dashed). 



obtained from the parameterization given in [19] (solid/black) to the threshold approxima- 
tion (I4.10p with the constant h q J set to zero. In the red/dot-dashed curve we multiply the 
correction (14.101) with the threshold approximation of the tree, in the blue/dashed curve we 
use the full tree cross section. It can be seen that the NLO corrections are peaked closer 
to threshold than the tree cross section, as might be expected. Multiplying the scaling 
function by the full tree leads to an improved agreement with the full NLO result. Overall 
the approximation of the full NLO result by the threshold expansion is quite good, in the 
sense that integrating the approximated partonic cross section captures the bulk part of 
the exact result. 



4.1.2 Production from gluon fusion 

For the gluon induced process the s-channel diagram is again P-wave suppressed, so the 
dominant contribution comes from the t- and w-channel diagrams and the quartic vertex. 
These diagrams give 



l(0) 



iA\l } {gg ->• q iR q jR ) = ig s Sij{T a \T a2 } a3a4 (e^g^e2) 



(4.11) 



with the transverse metric g^ v = g^ u — ^(n^n,^ + n^n u ). This form of the scattering 
amplitude at threshold is reproduced by the matrix element of the effective theory operator 
(no sum over flavours is implied) 



O (o) 



A 1 - A- 7/jt i/)' t 



(4.12) 
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Here it was used that the polarization vector corresponding to a free SCET gluon fields 
A 1 is given by e^_ = g ± ^ lJ e u . This can be seen from (I2.26P using a Fourier transformation 
after setting the strong coupling to zero, so that W c = 1 and the non-abelian terms in the 
field strength vanish. 

To introduce the colour basis, note that the gluon-gluon system admits the decomposi- 
tion 8® 8 = I©8 s ©8 a ffil0ffil0©27. The squark-antisquark pair is again either in a singlet 
or octet state, but this time the latter can be produced by two different initial colour-octet 
states. The colour basis obtained for this case from the prescription (13. 2p agrees with the 
one given in [TJ, 

(1) 1 r c 

° {a} V^Dl 3a4 

„( 2 ) _ jbaia^rpb 

< a > ' V2D A B F 



C {a} ~ V N£) ^0201-^1304' (4-13) 

where d abc are the usual symmetric invariant tensors of 577(3) and F% = if aiba2 . We also 



for the three colour channels are (including the non-relativistic normalization 2iriq) 



defined the coefficients D A = N% — 1 and B F = ° N = ^. Thus, the matching coefficients 



r (o,i) _ f^TT 2 Ka s , 



.9.9, \fJ-i v A r ^vj_ y Ml 



The diagonal elements of the hard function are given by 

#f = A^ f^Y -r^gf = ( ^V h^ (4.15) 

8m? \mq J AD\^ \ZD A m\) 

with the colour factors 

hf = 2C F , hf = 2D A B F , hf = . (4.16) 

Using Y2i hi = Cf(Nc — 2) gives for the total partonic cross section at threshold 

<x (0) I - ^ N "~ 2 U 17 ) 

gSHUkriJ 8 -*"* 16m? N C (N? - 1) ' { ' 

Adding the same expression for left and right squark labels exchanged again reproduces 
the result of [85]. In figure H] we compare the partonic cross section at threshold to the 
full partonic LO cross section [85], for a squark mass of 1 TeV at the LHC with 14 TeV 
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Figure 4: Partonic cross sections for gg — > qq multiplied with the gluon-gluon luminosity 
for rrig = 1 TeV at the 14 TeV LHC. Left panel: exact leading-order result (black/solid) 
and threshold approximation (blue/dashed). Right panel: exact NLO result (black/solid) 
and the approximation (I4.18J) with /% (/-0 = 0, using the threshold approximation of the 
tree (red/dot-dashed) and the full tree (blue/dashed). 



cms energy. As for the quark-antiquark induced subprocess we plot the product of the 

ninosity and the Jacobian ||. Again the threshold 
0.3 while the maximum of the integrand occurs 



partonic cross section with the parton luminosity and the Jacobian |5. Again the threshold 



approximation is adequate up to f3 
around j5 « 0.5. 

For the threshold expansion of the NLO cross section we obtain 
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+ hV>Qi) + 0(p). (4.18) 



The /3-dependent terms agree with [85J. Again, the hard corrections h gg are currently not 
available analytically. In the right plot of figure H] we compare the full NLO correction 
in the parameterization of pj5] (black/solid) to the threshold approximation of the NLO 
correction (14. 18ft with h gg = 0, using the threshold approximation of the tree (red/dot- 
dashed) or the full tree cross section (blue/dashed) as pre-factor. In this case, the threshold 
approximation gives a better agreement with the full NLO result than for the quark- 
antiquark initial state. Using the full tree further improves the agreement with the full NLO 
correction in the full /3-range. Once again the integrated threshold approximation provides 
a very good approximation to the exact NLO correction, especially when multiplying (14. 18ft 
by the full Born cross section. This can be compared to the Born cross section itself, for 
which the leading term in the /3-expansion is typically a less reliable approximation. 
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4.2 Results for squark-antisquark production at the LHC and 
Tevatron 

We are now ready to perform the combined soft and Coulomb NLL resummation for squark- 
antisquark production. We first discuss our choice for the soft, hard and Coulomb scales 
entering the resummed cross section, and then present numerical results for the LHC and 
the Tevatron. We also compare to the results of a soft gluon resummation in Mellin space 
presented in [20] . 

4.2.1 Scale choices 

Evaluating the NLL resummed cross section (13.531) requires a choice for the soft, hard and 
Coulomb scales. As can be seen from the Sudakov logarithm in the evolution equation of 
the hard function ( 13.331) . the natural scale of the hard corrections is of the order fih ~ 2M. 
Therefore we will choose the default value of the hard scale to be jlh = 2irtq. Analogously, 
the form of the approximate NLO corrections (14.101) and (14. 18j) implies that the scale elim- 
inating large logarithms in the soft corrections to the partonic cross section is given by 
fi s ~ 8E « 8M/3 2 . However, it was argued in [321 [33] that this choice leads to similar 
problems as the inversion of Mellin-space resummation formulae for the partonic cross sec- 
tion [TO], i.e. an ill-defined convolution with the parton luminosity. Here we follow [233 an d 
choose a soft scale such that it minimizes the relative fixed-order one-loop soft correction 
to the hadronic cross section. More precisely, we vary the scale in the PDFs and the soft 
correction and determine the value jl s that minimizes the relative soft corrections: 

d <s-^ f ati~ \^RP , ,eoft( rs >/M (a in \ 

0= w,^L v , dTV(T ^^w^r (4 - 19) 

Here the fixed-order NLO soft correction <5pJ soft is obtained from (14. 10J) and (14. 18j) . setting 
the Coulomb correction and the hard corrections hi to zero, and we divide by the leading- 
order hadronic cross section. This procedure results in values 

/} s = 123 - 455 GeV for m q = 200 - 2000 GeV (4.20) 

for the LHC with yfs = 14TeV. The fact that the soft scale is over half of the mass 
for light squarks, but considerably smaller for larger masses, indicates that the threshold 
resummation is more important for heavier squarks, as expected. We can approximately 
fit the numerical results for the soft scale by a function of the form proposed in 



~lhc w rn,(l - p) 

^ s V 3-8 + 149p ' [ ' 

with p = 4m~/s, which is accurate to better than 5% for rriq > 450 GeV. We find that for 
different cms energies \fs the ratio /2^ HC /rriq is to a good approximation a function of p 
only, that is, for a squark of mass 500 GeV at 7 TeV cms energy the ratio is approximately 
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Figure 5: Dependence of the NLL corrections at the 14 TeV LHC on the soft scale (left 
panel) and the hard scale (right panel), normalized to the corrections at the default scales 
fx a determined from (14.191) and jlh = 2rriq. Black (solid): rriq = 500 GeV, blue (dashed): 



rriq = 1 TeV, red (dot-dashed): m q 



2 TeV 



the same as for a 1 TeV squark at 14 TeV. In our numerical NLL results, however, we use 
the numerical values for the soft scale obtained directly by the condition (I4.19p . At the 
Tevatron the values for the soft scale can be fitted by the function 



-TeV 



rriq{l 



P) 



V1.7 + 81.3p' 
and representative values for the default soft scale are 
r (200GeV) = 84GeV , /tJ eV (400GeV) = 86GeV , 



(4.22) 



-TeV, 



/} s leV (600GeV) = 65GeV (4.23) 



Note that the soft scale begins to decrease above a certain value of p. 

Figure [5] shows the dependence of the NLL corrections to the hadronic cross section 
at the LHC with 14 TeV on the soft and hard scales, varied around the default values jl s 
and jlh = 2rriq. Here by "NLL corrections" we mean the convolution of the expression 
A<5"^ L (s, A) in (I3.55P with the parton luminosity, i.e. the corrections on top of the exact 
NLO cross section due to the higher-order soft-gluon and Coulomb terms (bound-state 
corrections are not included in these curves, and we set the cutoff A to the maximum 
possible value). In the left panel, the hard scale is fixed to the default value /i/i = jlh and 
we vary the soft scale in the interval 0.5/2 s . . . 2ji s . Analogously, in the right panel the soft 
scale is fixed to the default value and we vary the hard scale in the interval rriq . . . Arriq. 
The ambiguity due to the choice of the soft scale becomes smaller for larger masses, in 
agreement with the expectation that soft-gluon resummation is better justified for larger 
masses, where the contribution from the threshold region to the total cross section is more 
important. 

We now address the choice of the scale used for the strong coupling constant in the 
Coulomb Green function. The only scale dependence of the imaginary part of the leading- 
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Figure 6: Dependence of the NLL corrections at the 14 TeV LHC on the Coulomb scale, 
normalized to the corrections at the default scale fie — niax{2m g /3,/is} • Black (solid): 

2 TeV 



rriq = 500 GeV, blue (dashed): rriq = 1 TeV, red (dot-dashed): m 



order Coulomb Green function ( I3.27P enters through the coupling constant of the leading- 
order Coulomb potential, so the potential function is separately renormalization-group 
invariant at NLL, and the coupling constant can be evaluated at a scale fie different 
from the factorization scale. Since the necessary truncation of the perturbative series 
for the running coupling constant introduces a residual higher-order scale dependence, 
the scale \xc should be chosen such that higher-order corrections are minimized. Since 
the Coulomb corrections are related to the exchange of potential gluons with momentum 
transfer \k\ ~ Ma/A, a natural scale choice ixc is expected to be of the order 



fie ~ M\f\ = M(3 ~ Ma s 



(4.24) 



where we have used the NRQCD counting (3 ~ a s . Indeed, as mentioned in section [2.5.21 
the effect of the scale-dependent strength of the Coulomb potential can be incorporated 
using the choice \ic = 2rriqf3. The choice He ~ 0{Mj3) is in fact required to sum correctly 
all NLL terms, since He ~ M would miss terms such as ot 2 s /f5 x ln/3, which arise in part 
from the small virtuality of Coulomb gluons. However, very small values of /3 are integrated 
over in the convolution of the partonic cross section with the PDFs, and with the choice 
He oc (3 the strong coupling a s (nc) would hit the Landau pole. On the other hand, the 
relevant scale for bound-state effects is set by the inverse Bohr radius of the first HH' 
bound state, 1/rg = CFrriqa s /2 = /is/2. Therefore, as our default scale choice for the 
Coulomb Green function we use 



He = max{2rriqf3,HB} 



(4.25) 



where we solve the equation \xb = CFfnga s (fiB) iteratively. The Bohr scale was used in 
recent studies of threshold effects in top-pair or gluino production [241426] , while in the 
recent calculation [18] of squark-antisquark and gluino production the Coulomb scale was 
set equal to the factorization scale. The dependence of the NLL corrections at the LHC 
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with 14 TeV cms energy on the Coulomb scale (with all other scales fixed to their default 
values) is shown in Figure El where we vary the Coulomb scale in the interval 0.5/ic • • • 2fic- 
The scale dependence is of a similar magnitude as for the soft and hard scales, and again 
improves for larger masses. In the results given below, we estimate the uncertainty due 
to the scale choices by varying all scales /z/, /i s , \ih and \ic from one-half to twice their 
default value and add the uncertainties in quadrature. 

4.2.2 Results 

Having fixed the default scale values, we present numerical resultd_| for the combined NLL 
resummation of soft and Coulomb effects obtained by inserting the potential function ( I3.28P 
into the cross section (13.531) . We use the MSTW2008 PDF set |88j and the associated value 
of the strong coupling. The results for the fixed-order NLO cross section are obtained using 
the program Prospino, which is based on the calculation of [85]. In contrast to the default 
setting of Prospino, we include contributions from initial-state bottom quarks using b- 
quark PDFs. Unless stated otherwise, we use the squark-gluino mass ratio of m^ = 1.25rrig. 
As discussed in section 14.11 we consider the squarks to be mass-degenerate and do not 
include top squark final states. Our default value for the factorization scale is jlf = uiq. 

In table [1] we compare the Prospino NLO predictions obtained using the MSTW08 
and the CTEQ6.6M [HH] PDF set for the Tevatron and the LHC at 14 TeV. While we 
find good agreement for small squark masses, the differences between the two sets grow to 
over 10% for large masses. Since we expect very similar effects for the results including 
NLL resummation, we present our predictions using the MSTW PDFs below, but the 
uncertainties due to differences between current PDF sets at large x should be taken into 
account in interpreting our results. 

In figure [7] we show the NLL corrections to the partonic cross sections multiplied by the 
parton luminosities and the Jacobian |^. It can be seen that the dominant contributions 
arise from the threshold region f3 < 0.2, and that there is only a small difference between 
using the hard function obtained from the tree at threshold, as determined in section 14.11 
or from the full tree, as mentioned in section 13.21 Thus, in the results presented below we 
use the hard functions determined from the tree at threshold, and set the cutoff A in the 
matching formula ( 12. 17ft to the maximal value, so that the resummed correction is applied 
in the full phase space. Another plausible choice is to switch off the NLL corrections at 
the point where they turn negative [TJ]. In the example shown in figure [7] this occurs 
at /3 ~ 0.4, where the validity of the threshold approximation is not immediately clear. 
Since the integrated contributions from the region j3 > 0.4 are moderate and negative, our 
procedure gives a conservative estimate of the NLL corrections. We stress that this choice 
depends on a study of the behaviour of the resummed corrections at larger values of (3 and 
that other choices might be required for other processes. 



11 There are some differences to preliminary results presented in |50p51j due to a different choice of hard 
scale (/i/j = rriq) in the NLL results, and because the qg initiated subprocess was not included in the NLO 
cross section. 
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a(pp — > qq)(pb), y/s = 1.96 TeV 


m g -[GeV] 


NLOmstw 


NLOcteq 


200 


11.67{i 


11.633 


300 


6.687{ -^ x 10" 1 


6.93+t;^ x 10~ 2 


400 


4311SS* x 10- 


4.64^ x 10" 2 


500 


2.397a53 x l ^' s 


2.73+^ x lO" 3 





cr(pp — > qq)(pb), y/s = 14 TeV 


rrig[GeV] 


NLOmstw 


NLOcteq 


500 


14.435 


13.73- 


1000 


2.587^? x 10- 1 


2.497a29 x 10- 1 


1500 


1.357^ x 10- 2 


1.383$ x lO" 2 


2000 


9.977{;l4 x 10~ 4 


1.137^6 x 10- y 



Table 1: NLO results for the squark-antisquark production cross section at the Teva- 
tron and the LHC at y/s = 14 TeV for m g = 1.25m, obtained with Prospino, using 
MSTW08NLO and CTEQ6.6M PDFs. The error estimate is obtained by varying the 
factorization scale in the interval rriq/2 < [if < 2m g - 
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Figure 7: NLL corrections to the partonic cross sections, multiplied with the parton lu- 
minosity, for qtfj — >■ %qi (left) and gg —¥ %qi at the 14 TeV LHC for m, = 1 TeV and 
rrig = 1.25 TeV, using the threshold approximation of the tree (red, dot-dashed) and the 
full tree (blue/dashed). 
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Figure 8: NLL corrections to the squark-antisquark cross section at the LHC with y/s = 
14 TeV (above) and the Tevatron (below) relative to the NLO cross sections in various 
approximations. For explanations see text. 



In figure [8] we plot the i^-factor relative to the NLO corrections, 



AK, 



NLL 



C"NLL ~ °"NLO 
C"NLO 



(4.26) 



for the y/s = 14 TeV LHC and for the Tevatron, as a function of the squark mass, and 
consider various approximations to our NLL result to study the effect of the Coulomb and 
soft corrections separately. The different curves in the figure are defined as follows: 

NLL (black, solid): Results for the full NLL soft and Coulomb resummation using the 
resummed cross section ( 13.531) in the matching formula ( 13.551) . including bound-state 
effects obtained from the potential function below threshold ( I3.30p . 

NLL noB s (green, dashed): Bound-state effects are omitted from the full NLL result. 
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NLL s _|_h (blue, dashed): NLL resummation without Coulomb corrections, i.e. the cross 

2 / 

section ( I3.53[) with the trivial potential function J^(E) = -^\ N^z- 

NLL s _|_h+C (red, dot-dashed): The pure Coulomb corrections, obtained by using the 
trivial soft function W^ ' = 5(co) and the resummed potential function (I3.28J) without 
bound state effects, are added to the NLL^ approximation. 

It can be seen that the soft-Coulomb interference effects, given by the difference of the 
NLL noB s and NLL s+h +C curves, are sizable, in particular for large squark masses. Most 
of this effect is due to the interference of the first Coulomb correction with the resummed 
soft corrections; for instance at the 14 TeV LHC the relative correction at rriq = 2 TeV in 
this approximation is AKnll = 9.0% as opposed to AT^nll = 9.7% in the full NLL soft 
plus Coulomb resummation. The NLL s+ h corrections are of a similar magnitude as the soft 
NLL corrections obtained in [T7llT8~t l2"0"] from resummation in Mellin space. A quantitative 
comparison to these previous results will be performed below. The Coulomb corrections 
are larger than those observed in [18]. This can be traced to the choice of the Coulomb 
scale, which was taken to be \iq — A 4 / i n that reference, while we use the smaller scale (I4.25J) 
required to sum all NLL terms, as discussed above. 

In tables [2] and [3] we compare the NLL resummed results pertaining to pp collisions at 
the Tevatron cms energy -y/i = 1.96 TeV and pp collisions at LHC with -y/i = 7, 10 and 
14 TeV to the LO and NLO results obtained from Prospino. For the LO predictions we 
have used the MSTW2008LO set of PDFs, for the NLO and resummed predictions the 
MSTW2008NLO set. One observes a sizable reduction in the scale dependence for larger 
squark masses. This improvement can also be seen in figure [9j In the left-hand plots the 
scale dependence of the NLO cross section, the full NLL and the NLL s+ h correction are 
plotted as a function of the squark mass for the Tevatron and the 14 TeV LHC. For the 
NLL results, the soft, hard, Coulomb scales and the factorization scale are varied and the 
uncertainties added in quadrature. Clearly, a significant reduction of the scale dependence 
requires the inclusion of soft-Coulomb interference. In the right-hand plots in figure [9] 
we plot the LO, NLO and NLL cross section as a function of the factorization scale for 
rriq = 1 TeV at the 14 TeV LHC and for rriq = 400 GeV at the Tevatron. One observes a 
stabilization of the result with respect to variations of the factorization scale in going from 
the NLO to the NLL approximation if the factorization scale is varied in the usual interval 
0.5m,j < fif < 2rriq, even after taking the variation of the soft, hard and Coulomb scales 
(shown as the green band) into account. 

In table H] we compare to the NLL result of [20] obtained from threshold resummation 
in Mellin space. In that reference no Coulomb resummation was considered. In order to 
compare to their result we use the approximation NLL s+ h defined above, but in addition 
set the hard scale to \Xh, — (J*f- As can be seen from the comparison, this approximation, 
denoted by NLL S , is in good agreement with the result from [20] at the 14 TeV LHC, 
while the full NLL corrections including Coulomb effects are considerably larger at higher 
masses. At the Tevatron the NLL S approximation results in smaller corrections than the 
ones obtained in the Mellin space approach, while again the Coulomb corrections and 
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a(pp — ¥ qq)(pb), y/s = 1.96 TeV 




m g -[GeV] 


LO 


NLO 
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A^ NLL 


200 


y.oy_ 2 .76 


ii-6ii:S 
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500 


1.87™ x 10- 
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cr(j)p — > qq)(pb), y/s = 7 TeV 




m g ~[GeV] 


LO 


NLO 
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Ai^NLL 
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<r(pp — > qq)(pb), y/s = 10 TeV 




m g ~[GeV] 


LO 


NLO 


NLL 


AiT NLL 
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1000 
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5.32™ x IO- 2 
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8.9% 
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1.02™ x IO- 2 


1.18™ x IO- 2 


1.30-;- x 10- 


10.6% 


1400 


2.48+^ x IO" 3 


2.79™ x 10- 


3.15^ x IO" 3 


12.7% 



Table 2: Results for the squark-antisquark production cross section at the Tevatron and 
the LHC at \fs = 7 TeV and 10 TeV for m g = 1.25rrig at leading order, next-to- leading 
order (NLO) and NLO plus resummed soft and Coulomb corrections (NLL). The error 
estimates of the LO and NLO results are obtained by varying the factorization scale in the 
interval rriq/2 < fif < 2m,q, the error estimate of the NLL result is obtained by varying 
Hi G {fif, fih, Us, He} i n the interval jli/2 < //» < 2/Zj around the default values and adding 
the variations of the cross section in quadrature. 
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cr(jpp -> gg)(pb; 


, VS = 14 TeV 
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Table 3: Results for the leading order, fixed-order NLO and resummed soft + Coulomb 
corrections (NLL) to squark-antisquark production at the LHC at y/s = 14 TeV, setup as 
in table El 





a{pp — > qq)(pb), a/s = 14 TeV 


m,-[GeV] 


NLO 


NLL M eiiin (ref. [20J) 
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1.61 x 10 1 (1.2%) 


1.62 x IO 1 (1.3%) 


1.67 x IO 1 (4.2%) 


1000 


2.89 x 10- 


2.93 x IO" 1 (1.7%) 
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a(pp —¥ qq)(pb), yfs = 1.96 TeV 


m g ~[GeV] 


NLO 


NLL M ciiin (ref. [20j) 


NLL S 


NLL 


200 
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9.79 x 10- 
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Table 4: Comparison to the results of [20] for the LHC at 14 TeV (above) and the Tevatron 

(below) for m~ g 

result. 



rrig. In brackets we quote the corrections A/^nll relative to the NLO 
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Figure 9: Left: Scale dependence of the NLO, NLL and NLL s+ h approximations at the 
LHC with y/s = 14 TeV (top) and the Tevatron with \/s = 1.96 TeV (bottom). The 
corrections are normalized to the value at the default scales. Right: Dependence of the 
LO, NLO and NLL result on the factorization scale, for rrtq = 1 TeV at the 14 TeV LHC 
(top) and for rriq = 400 GeV at the Tevatron (bottom). The green band indicates the 
uncertainty on the NLL result from the choice of the soft, hard and Coulomb scales. 



mixed soft/gluon corrections are sizable (note that, due to statistical fluctuations, some of 
the NLO results from Prospino differ from the numbers quoted in [20] in the last digit). 

5 Conclusion 

The production of pairs of massive coloured particles is an integral part of hadron collider 
physics programmes, since such particles, the top quark for certain, may be produced in 
large numbers. Precise calculations of the production cross sections are hampered by po- 
tentially large quantum corrections originating from the suppression of soft gluon emission 
and the exchange of Coulomb gluons near the partonic production threshold, especially for 
large masses of the produced particles. The techniques for summing soft gluon effects to 
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all orders in perturbation theory have been available for some time, and have been used 
extensively to improve Drell-Yan production and related processes. Likewise, the summa- 
tion of Coulomb effects is well established in the context of calculations of top quark pair 
production near threshold in e + e~ collisions. But the interplay of both effects as relevant 
to hadronic pair production had not been studied to all orders up to now. 

Joint soft-gluon and Coulomb resummation is a non-trivial issue, since the energy of 
soft gluons is of the same order as the kinetic energy of the heavy particles produced. 
Soft-gluon lines may therefore connect to the heavy-particle propagators in between the 
Coulomb ladders, as well as to the Coulomb gluons itself, impeding the standard factoriza- 
tion arguments that assume either no couplings to the final state, or to energetic particles. 
In this paper we employed a combination of soft-collinear and non-relativistic effective field 
theory techniques to show that, despite these apparent complications, soft gluons can be 
decoupled from the non-relativistic final state. Their effect is then contained in a more 
complicated soft function including two time-like Wilson lines for the final state particles, 
and in a convolution in energy of the soft function with the final state Coulomb function, 
since the latter depends on the energy but not the momentum taken away by the radiated 
gluons. Our main result is therefore a new factorization formula for the partonic cross sec- 
tion of heavy-particle pair production near threshold in hadron collisions given by (I2.63p . 
and the corresponding equation ( I2.93P for the invariant mass distribution, together with 
the resummation formula f l3.48|) for the partonic cross section. 

The result is naturally formulated in momentum space. When expressed in Mellin 
moment space, the factorization is multiplicative in each separate colour channel, and given 
by the product of a hard, soft and Coulomb (or potential) function. This justifies previous 
treatments of Coulomb effects at the one-loop order and extends them to all orders. The 
resummation of Coulomb effects to all orders is in fact formally required even at the leading- 
logarithmic level, though numerically the two- and higher-loop Coulomb correction is rather 
small. Starting from the NNLL level the combined soft-gluon and Coulomb resummation 
might have involved subleading soft gluon couplings from the effective Lagrangians (beyond 
the standard eikonal approximation in the collinear sector), making the resummation much 
more difficult. We applied symmetry considerations to show that these couplings do not 
contribute at NNLL. A new effect at this order is, however, that the potential function 
acquires factorization-scale dependence that cancels with the hard function, resulting in 
additional threshold logarithms on top of the ones from soft gluon emission. Eq. (I2.63P is 
thus valid at NNLL, but beyond that order further, more complicated soft functions are 
most likely required. 

Factorization and resummation are largely model-independent. The details of the hard 
production mechanism enter only in the momentum-independent short-distance coeffi- 
cients. We exemplified our theoretical approach by considering squark-antiquark produc- 
tion in pp and pp collisions at LHC and Tevatron energies with NLL accuracy. We find that 
the effect of resummation on top of the fixed-order NLO cross section is around (12-27)% 
for squark masses that give cross sections at the femtobarn level, and increases with mass 
as expected. The scale uncertainty is considerably reduced after resummation suggesting 
a reduction of the remaining theoretical error from which experimental limits or determi- 
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nations of squark masses would benefit. We find resummation effects larger than previous 
estimates in the literature, since we normalize the Coulomb function at the scale Mf3 to 
sum all NLL terms, and include the interference terms between soft gluon and Coulomb 
exchange, as well as squark-antisquark bound-state production in the total cross section. 
A more extensive phenomenological analysis of sparticle production covering all final states 
at the NNLL level is therefore of interest. 
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A Relation of collinear matrix elements to PDFs 

In this appendix we show that collinear matrix elements of the form (p\(ft c (x)<j) c (0)\p) can 
be expressed in terms of the parton distribution functions (PDFs) and a polarization sum 
of the external wavef unctions, as stated in (12.521) . We recall the operator definitions of the 
PDFs of quarks and gluons in full QCD |90j : 

f q/N (x x ,ti) = ^ J dte-^ m ^{NMtn)^W n (tn,0)^(0)\N) , (A.l) 

f 9/N (x 1 ,n) = X _ - — J dte- ix ^(N\{n,F^{tn)W n {ntMn P Fv P m\N) , (A.2) 

where ip is the quark field in QCD, F^ v the full gluon field strength and an average over 
the physical polarization states is understood. W n (x + ,0) denotes a Wilson line with the 
full gluon field n ■ A extending from to x + , 



W n (x + ,0)=Pexp 



ig s J dtn ■ A(nt) 
o 



(A.3) 



The PDFs in SCET have been introduced in [52] and are collected in [HI] in the SCET 
conventions used here. The purpose of the present discussion is to fix our conventions and 
identify all the numerical prefactors. 

In order to relate the collinear matrix elements to the PDFs, we decompose the collinear 
and anticollinear matrix elements into a basis of spin structures T l , so that (suppressing 
colour indices) 

(p\4>U<f>c;a\p) = J2 f U (Pl^TVclp) (A.4) 
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with tr[PP] = 1. We now discuss the cases of quarks and gluons separately. For n-collinear 
quarks the fields are given by C = W^ c . The identities r/£ c = and {rjifji '■/ '4)£ c = £ c 
satisfied by the collinear spinors imply [S2] that the independent spin structures for the 
decomposition of a matrix element (p|£ c r£c|p) can be taken as P G {^,^7 5 ,^7^}, up to 
normalization. Since we sum over the initial state parton spin, only the ff, term contributes. 
In colour space the product of £' and £ decomposes into a singlet and an octet component 
but the sum over the initial state parton colours projects on the colour singlet. We therefore 
obtain for the collinear quark matrix element 



<.(°)t/^(o) 



{p{Pi)\€^)€m\p(Pi))\ 



avg. 



= ±-6* (f/) {p|[e ) Wc (0) ](^)|[W c (O)t e ) ](0)b}| avg . 

= KpWhj I dx 1 e lx ^ p ^f q/p (x 1 ^). (A.5) 

Jo 

In the last step we have identified the quark PDF defined in terms of the SCET fields 

/(,/pK/*) = ^Jdte-^^ (p|[£l o) W c (o) ](^)|[Wi o )t£(°)](0)b)| avg . , (A.6) 

which is consistent with the definition f |2.52j) in the main text, since for quarks the spin- 
dependent normalization factor, identified as a polarization sum of collinear quark spinors 
with momentum P£ = (n ■ Pi)/2n M , 

*i(*> - ^ (f ') = ^ £ *(*)<(«) • (A.7) 

satisfies N q {x\P{) = xiN q (Pi). (The factor 7 arises because of the definition of the 
conjugate quark and antiquark fields as discussed below (12.251) .) Eq. (1A.6I) is of the same 
form as the PDF in full QCD (lA.ip for external quark states, apart from the fact that 
the collinear quarks and Wilson lines appear instead of the full QCD fields. Since we are 
considering external massless on-shell partons with vanishing transverse momentum, the 
collinear region is the only one contributing to the PDFs. As the SCET Lagrangian for a 
single collinear direction is equivalent to full QCD [59], one can replace the collinear fields 
in (JA.6P by the full QCD fields and identify the PDF with flA.ip . For a matrix element 
with collinear antiquarks we obtain an analogous expression with the antiquark PDF 

/*/,(*!, ^ = ^J dt e " M ' Pl (P\ tr { I [^ c (0)t ^ 0) ] («0 & 0) W^} (0) } |p> I spin avg . (A.8) 
and the normalization factor 
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For the matrix element of n-collinear gluons, the only available transverse symmetric 
tensor is g^ u so that [B"2"|l*JTj 

(p\^c;k^ + )Aj. iu (0)\p)\ avg . = 2 < N2 _ jt gf LU fe(p|^ M (z+)^ > (0)|p)| avg . 

= N^(P l )5 kl ^/ 9/p (x 1;/ u)e^^) (A.10) 

Jo x i 

with the gluon PDF in SCET 

f a/p (x ul i) = - ^f^ J tie-**** (p\A^(tn)AJjO)\p)U, (All) 

which is consistent with the definition f |2.52j) in the main text. The normalization factor is 

wx) - 5^ - m=T)Z<™*™ ■ (A - 12) 

For the last identity we use the polarization sum of the gluon polarization vectors 

2^ s( p i) e - ( p i) = -^ + — ^— — = -g*» + — — - — ( A - 13 ) 

r \ • q u ■ q 

and choose the arbitrary light-like vector q proportional to n. 

The equivalence of the SCET and QCD definitions of the gluon PDF is established 
most easily in light-cone gauge n • A = 0. In this gauge the SCET gluon operators (I2.26P 
reduce to the collinear gluon fields: Ajr = A^r, while the QCD definition (I A. 21) becomes: 

/,/p(si, Ai) = ^—^t [dte-^ tn - p i (p\(n ■ dA»(tn)(n ■ &4„(0))b>j avg . . (A.14) 

27rxi{n ■ F\) J 

Upon integration by parts and Fourier transformation one sees that this agrees with the 
SCET definition. 

B Coulomb potential for gluinos and squarks 

Here we collect the coefficients D Ra of the Coulomb potential relevant for the production 
of coloured SUSY particles defined by the relation (13.201) . This relation can be solved for 
arbitrary representations in terms of the quadratic Casimir operators of the representations 
of the two heavy particles and the final state system 

D Ra = ~(C Ra -C R -C R >), (B.l) 
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as can be deduced from the generator of the product representation r £( R ® R ) = ^(R) ® 
l(R') _|_ ]_(-#) (g) ^{R') anc [ by projecting the identity 

T (fl)6 T (H')6 = I I TmR ')b T{m Ri)b _ ( Cr + CR ,)lW ® l(«')l (B.2) 

on the irreducible representation. The explicit results for gluino pair production are long 
known [H2] while all the remaining cases have recently been collected in [93J. We provide 
these results here for completeness. Note that our different sign convention implies that 
negative values of Du a correspond to an attractive Coulomb force, positive values to a 
repulsive one. 

S quark- antis quark production. The coefficients for squark-antisquark production have been 
quoted already in (13. 23ft : 



n - r - N c~ l - 4 n Ca r l l (W 

U 1 --C F - ___-_ d 8 _ -c F ____-. (B.S) 



^^ = -- D 8 = ^ 

2N C 3 ' 2 2N C 6 



Squark-squark production. In squark-squark production qq — > qq the final state system is 
either in the 3 or 6 representation and the coefficients of the Coulomb potential are 

D5 =-K i+ i)=-5- D6 =K i_ i) = 5' (R4) 

Gluino-squark production. For gluino-squark production qg — > qg the final-state represen- 
tations appear in the decomposition 3©8 = 3 + 6 + 15. The coefficients read 

d* = -Y = ~1> ^ = 4' Dl " = + \- (R5) 

Gluino-pair production. For gluino pair production the final state representations appear 
in the decomposition 8<g}8 = l©8 s ©8 a ffil0ffil0ffi27 with the coefficients of the Coulomb 
potential: 

D x = -3, D S3 =D 8a = ~, (B.6) 

D 10 = 0, D 27 = l. (B.7) 

Note that only the singlet and octet states can be produced from a quark-antiquark initial 
state while all the states appear in gluon fusion. 

C Evolution functions 

For convenience we quote here the explicit expressions of the resummation functions S, 
a( and a? given in ( I3.42p up to NLL. The expressions for the corresponding functions for 
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the Drell-Yan process up to N 3 LL order have been given in [33]. Up to NLL the relevant 
expressions read 



S(ji h ,n, 



r (0),r , r (0)r' - 

J- CUSD 1" A c 




V a s (fJ> s ) a s (fj, h )J 




a s (fi s ) , . a s (fi s )\ X 2 a,(|J,) 

-; r + in 7 r -\ — in -. r 


, (C.l) 


•) 


(C.2) 




(C.3) 



2/3 a 8 (fi h ) 

and an analogous expression for a^ ,r . Here we have introduced the perturbative expansions 
of the beta function and the anomalous dimensions: 

' Os_ \ n+1 

Air, 



/3(a s 
7(«« 



-2a s ^ /3„ Q 



n=0 



£~ (B) 



,7 V 



n=0 
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Air 



(C.4) 
(C.5) 



with 



0o 






0i 



~r 2 - 



10 



C471/ — ICfUj . 



(C.6) 
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T L/i 3 

Explicit expressions for the anomalous dimensions required in the present work can be 
found, e.g. in 



D Fixed-order expansions 



The expression for the resummed cross section (I3.48P can be expanded to a fixed order 
a™ in the strong coupling, providing an approximation to the full 0(a™) QCD calculation. 
According to the counting (12.71) . the expansion of the NLL resummed cross section to order 
a s compared to the leading order cross section is accurate up to terms of the order a s x 0°. 
This approximation to the NLO cross section is obtained from (13.531) by expanding the 
evolution function [/, to 0(a s ), performing the convolution of the leading order Coulomb- 
Green function for arbitrary r], identifying i] = 2a r a (fj, s , /x/) and expanding up to 0(a s ). 
Inserting the explicit results for the one-loop cusp and soft anomalous dimensions, the 
result reads 



*NLL(1) 

J pp'(i) 



2n 2 D R 

Ha 



B 



"T y^r t Ly 



+ ln 2 
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In 2 
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'(SH©) 


4. 
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M 2 Y 
*f - 


- 4 (C^ + 4 (CV + CVO) In f— 
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+ 2 ( 7 ^(°) + 7 ^'(°) + 2C Ra - 2ft) In f^\ + (9(1) . (D.l) 

We have introduced the expansion of the cross section in the strong coupling 

<V(p> Hf) - <v < i + -^p/^w + l^^^J V« + u t a J [ • (DJ) 

The term involving 2/9 in the last line of (ID.lj) appears, since we express o pp , oc cv^ 
in terms of the strong coupling at scale [if, while in the factorization formula the hard 
function involves the scale jMh- The expansion (ID.lj) is used in (I3.55P to match the NLL 
resummed prediction for squark-antisquark production to the fixed-order NLO calculation. 
At the order considered here, one can further approximate E = MP 2 + 0(f3 A ) to write the 
result in a more familiar form. 

The NLL approximation to the NLO cross section (ID.lj) can be improved to include 
in addition all terms of order a s x /3° by inserting the Laplace transform of the one-loop 
soft function (13. 9p as well as the hard function in (13 .48 p . Performing the derivatives with 
respect to 77 and expanding to 0(a s ) afterwards we obtain for the NLO corrections to the 
cross section in the colour-channel % |45j : 



a„ = -^ + 4(C r + C,) 



Vw/ 24 



- 4 (Cr, + 4 (C r 4- CV)) In (— ") + 12C H „ + fcfV/) 4" 0(/3) ■ (D.3) 

We note that the dependence on the soft scale has canceled between the one-loop soft func- 
tion and the expansion of the evolution function. Similarly the evolution equation ( I3.33P 
implies that the dependence on the hard scale cancels between h\ (///J = H\ \jjLh)/H\ 
and the evolution function. This can also be checked explicitly for the case of top-pair 
production where the one-loop hard functions can be obtained from [80]. In contrast, the 
NLL result (jD.ip contains residual dependence on the soft and hard scale that is formally 
of higher order if the scales are chosen of the order /j, s ~ E ~ M(3 2 and \Xh ~ 1M. We 
observe that the In 8 constant terms in the NLO result (1D.3|) can be reproduced exactly 
with \x a = 8E. 

Expanding the NNLL resummed cross section in the same way results in an approx- 
imation to the NNLO cross sections fp P that is accurate up to terms a 2 /3°. The results 
for the production of heavy particles of arbitrary spin and colour and the application to 
top-pair production have already been presented in [45] and will not be repeated here. 

References 

[1] G. Sterman, Nucl. Phys. B281 (1987) 310. 



65 



[2] 
[3] 
[4] 

[5] 

[6] 

[7] 



S. Catani and L. Trentadue, Nucl. Phys. B327 (1989) 323. 



E. Laenen, J. Smith, and W. L. van Neerven, Nucl. Phys. B369 (1992) 543-599. 



E. L. Berger and H. Contopanagos, Phys. Rev. D54 (1996) 3085-3113 
|arXiv : hep-ph /9603326J 

S. Catani, M. L. Mangano, P. Nason, and L. Trentadue, 
\Nucl. Phys. B478 (1996) 273-310[ |arXiv : hep-ph /960435l| 



N. Kidonakis and G. Sterman, Phys. Lett. B387 (1996) 867-874 



N. Kidonakis and G. Sterman, \Nucl. Phys. B505 (1997) 321-348 
arXiv : hep-ph/9705234~] 



N. Kidonakis, G. Oderda, and G. Sterman, Nucl. Phys. B531 (1998) 365-402 
arXiv : hep-ph/ 980324T| 



[9] 
[10] 

[11] 
[12] 

[13] 

[14] 
[15] 

[i6; 

[17 
[18 
[19 



R. Bonciani, S. Catani, M. L. Mangano, and P. Nason, 



Phys. Lett. B575 (2003) 268-278, arXiv:hep-ph/0307035 



S. Catani, M. L. Mangano, P. Nason, and L. Trentadue, 



Phys. Lett. B378 (1996) 329-336, arXiv:hep-ph/9602208 



R. Bonciani, S. Catani, M. L. Mangano, and P. Nason, 
\Nucl. Phys. B529 (1998) 424-450HarXiv:hep-ph/9801375| 



N. Kidonakis, E. Laenen, S. Moch, and R. Vogt, \Phys. Rev. D64 (2001) 114001 , 
arXiv : hep -ph/0105041| 

M. Cacciari, S. Frixione, M. L. Mangano, P. Nason, and G. Ridolfi, JEEP 04 (2004) 
068, larXiv : hep-ph/0303085J 

S. Moch and P. Uwer, \Phys. Rev. D78 (2008) 0"34 003, a rXiv: 0804. 14761 Ihep-phTl 

M. Cacciari, S. Frixione, M. L. Mangano, P. Nason, and G. Ridolfi, 



JEEP 09 (2008) 127, arXiv: 0804. 2800 [hep-ph] 



N. Kidonakis and R. Vogt, \Phys. Rev. D7 8 (2008) 074005, 
larXiv: 0805 .38441 I hep-phi I 



A. Kulesza and L. Motyka Jj%g, Rev. Lett. 102 (2009) 111802 
larXiv: 0807 .24051 I hep-ph 1 1 



A. Kulesza and L. Motyka, Phys. Rev. D80 (2009) 095004, 
larXiv: 0905 .47491 I hep-phi 1 

U. Langenfeld and S.-O. Moch, Phys. Lett. B675 (2009) 210-221, 
arXiv : 0901 . 0802 [hep-ph] . 



66 



[20] 
[21] 

[22] 
[23] 
[24] 

[25] 

[26] 
[27] 

[28] 

[29] 
[30] 
[31] 
[32] 
[33] 

[34] 

[35] 

[36] 

[37] 

[38] 



W. Beenakker et al., J HEP 12 (2009) 041[ jrXiT: 0909 .4418 [hep-ph] , 



A. Idilbi, C. Kim, and T. Mehen, Ph ys. Rev. D79 (2009) 11 4016 
larXiv: 0903 .36681 Ihep^phTl 



A. Idilbi, C. Kim, and T. Mehen, arXiv: 1007.0865 [hep-ph] . 

V. S. Fadin, V. A. Khoze, and T. Sjostrand, \Z. Phys. C48 (1990) 613^622 



K. Hagiwara, Y. Sumino, and H. Yokoya, Ph ys. Lett. B666 ( 2008) 71-76 
larXiv: 0804 .10141 I hep-phi I 

Y. Kiyo, J. H. Kiihn, S. Moch, M. Steinhauser, and P. Uwer, 



Eur. Phys. J. C60 (2009) 375-386, arXiv: 0812. 0919 [hep-ph] 



K. Hagiwara and H. Yokoya, JEEP 10 (2009) 049, arXiv: 0909. 3204 [hep-ph] 



H. Contopanagos, E. Laenen, and G. Sterman, \Nucl. Phys. B484 (1997) 303-330 , 
|arXiv : hep-ph/96043 13] 

G. P. Korchemsky and G. Marchesini, Nucl. Phys. B406 (1993) 225-258, 



hep-ph/9210281, 



G. P. Korchemsky and G. Marchesini, Phys. Lett. B313 (1993) 433-440. 
A. V. Manohar, Phys. Rev. D68 (2003) 114019, hep-ph/0309176. 



T. Becher and M. Neubert, Phys. Rev. Lett. 97 (2006) 082001, hep-ph/0605050 



T. Becher, M. Neubert, and B. D. Pecjak, JEEP 01 (2007) 076, hep-ph/06 07228 



T. Becher, M . Neuber t, and G . Xu, JEEP 07 (2008) 030 
larXiv: 0710 .06801 Ihep^phTl 



V. Ahrens, T. Becher, M. Neubert, and L. L. Yang, \Phys. Rev. D79 (2009) 033013 
arXiv : 080 8.30081 I hep-ph ll 

V. Ahrens, T. Becher, M. Neubert, and L. L. Yang, 

\Eur. Phys. JTC 62 (20 09) 333-353] larXiv: 0809. 42831 Ihep^phTl 

V. Ahrens, A. Ferroglia, M. Neubert, B. D. Pecjak, and L. L. Yang, 
larXiv: 1003.58271 Ihep^phTl 



M. Beneke, P. Falgari, and C. Schwinn, \Nucl. Phys. B828 (2 010) 69-101 
larXiv: 0907. 14431 I hep-phi I 



T. Becher and M. Neubert, Phys. Rev. D79 (2009) 125004 



arXiv: 0904. 1021 fhe p-phl 



67 



[39] N. Kidonakis, \Phys. Rev. Lett. 102 (2009) 232003, arXiv: 0903 .2561 [hep-ph] . 



[40] G. P. Korchemsky and A. V. Radyushkin, Phys. Lett. B279 (1992) 359-366, 



arXiv : hep-ph/9203222 



[41] M. Czakon, A. Mitov, and G. Sterman, \Phys. Rev. D 80 (2009) 074017, 
larXiv : 0907 . 17901 fhep-phl I 

[42] A. Ferroglia, M. Neubert, B. D. Pecjak, and L. L. Yang, 

\Phys. Rev. Lett. 103 (2009) 20160T| larXiv : 0907 . 479 1 1 I hep^nTl 



[43] A. Ferroglia, M. Neubert, B. D. Pecjak, and L. L. Yang, JEEP 11 (2009) 062, 
larXiv: 0908 .36761 fhep-phl I 

[44] A. Mitov, G. Sterman, and I. Sung, larXiv: 1005~464 6 fhep ^phTl 

[45] M. Beneke, M. Czakon, P. Falgari, A. Mitov, and C. Schwinn, 



Phys. Lett. B690 (2010) 483-490, arXiv: 0911 .5166 [hep-ph] , 



[46] M. Beneke, A. P. Chapovsky, A. Signer, and G. Zanderighi, 



Phys. Rev. Lett. 93 (2004) 011602, arXiv:hep-ph/0312331 



[47] M. Beneke, A. P. Chapovsky, A. Signer, and G. Zanderighi, 
\Nucl. Phys. B686 (2004) 205-247| |arXiv : hep-ph/0401 002 , 

[48] M. Beneke, N. Kauer, A. Signer, and G. Zanderighi, Nucl. Phys. Proc. Suppl. 152 
(2006) 162-167, fjrXiv : hep-ph/041 1008| 

[49] M. Beneke, P. Falgari, C. Schwinn, A. Signer, and G. Zanderighi, 
\Nucl. Phys. B792 (2008) 89-135] larXiv: 0707. 07731 I hep-phi I 

[50] M. Beneke, P. Falgari, and C. Schwinn, PoS (EPS-HEP 2009) (2009) 319, 
larXiv: 0909. 34881 Ihep-phTl 

[51] M. Beneke, P. Falgari, and C. Schwinn, PoS (RADCOR2009) (2010) 012, 
larXiv: 1001.46271 fhep-phl 1 

[52] H. Baer, V. Barger, A. Lessa, and X. Tata, | JEEP 09 (2009) 063[ 
larXiv: 0907. 19221 fhep-phl 1 



[53] H. Baer, V. Barger, A. Lessa, and X. Tata, JEEP 06 (2010) 102 
larXiv: 1004.35941 fhep-phl 1 

[54] DO collaboration Collaboration, V. M. Abazov et al., 

\Phys. Lett. B660 (20 08) 449-4~57[ larXiv :0712~3 805 f hep-exl 1 

[55] CDF collaboration Collaboration, T. Aaltonen et al, 

\Phys. Rev. Lett. 102 (2009) 121801| larXiv : 081 1.25121 I hep-ex II 

68 



[56] 
[57] 
[58] 
[59] 
[60] 
[61] 
[62] 
[63] 
[64] 
[65] 



V. Ahrens, A. Ferroglia, M. Neubert, B. D. Pecjak, and L. L. Yang, 
\Phys. Lett. B687 (2010) 331-337J larXiv : 0912 .33751 Ihep-phll 

C. W. Bauer, S. Fleming, D. Pirjol, and I. W. Stewart, Phys. Rev. D63 (2001) 



114020, hep-ph/0011336 



C. W. Bauer, D. Pirjol, and I. W. Stewart, Phys. Rev. D65 (2002) 054022, 



|hep-ph/010 9045 



M. Beneke, A. P. Chapovsky, M. Diehl, and T. Feldmann, Nucl. Phys. B643 (2002) 



431-476, hep-ph/0206152 



M. Beneke and T. Feldmann, [Phys. Lett. B553 (200 3) 267-276 , 
arXiv : hep-ph/021 1358 , 



S. Actis, M. Beneke, P. Falgari, and C. Schwinn, Nucl. Phys. B807 (2009) 1-32, 
larXiv: 0807 .01021 Ihep-phll 

C. W. Bauer, S. Fleming, D. Pirjol, I. Z. Rothstein, and I. W. Stewart, Phys. Rev. 
D66 (2002) 014017, |hep-ph/02 02088. 

M. Beneke, P. Falgari, and C. Schwinn, PoS (RADCOR2009) (2010) 011, 
larXiv: 1001 .46211 rhep^phTl 

P. Falgari, W-pair production near threshold in unstable-particle effective theory. 
PhD thesis, RWTH Aachen University, 2008. 

T. Becher, R. J. Hill, B. O. Lange, and M. Neubert, Phys. Rev. D69 (2004) 034013, 



hep-ph/0309227. 



M. Beneke and V. A. Smirnov, \Nucl. Phys. B522 (1998) 321-344 
|arXiv :hep-ph/971139T 



[67] 

[68] 
[69] 

[70] 
[71] 



A. Pineda and J. Soto, \Nucl. Phys. Proc. Suppl. 64 (1998) 428-432, 
arXiv : hep-ph/9707481 



M. Beneke, arXiv:hep-ph/9806429 



N. Brambilla, A. Pineda, J. Soto, and A. Vairo, \Nucl. Phys. B566 (2 000) 275 
arXiv : hep-ph/9907240 . 



M. Beneke, arXiv:hep-ph/99 11490 Proceedings of the 8th International 
Symposium on Heavy Flavor Physics (Heavy Flavors 8), Southampton, England, 
25-29 Jul 1999. 



P. Labelle, Phys. Rev. D58 (1998) 093013, arXiv:hep-ph/9608491 



69 



[72; 

[73] 

[74] 

[75] 
[76] 

[77] 

[71 

[79] 

[80" 



B. Grinstein and I. Z. Rothstein, \Phys. Rev. D57 (1998) 78-82 
|arXiv : hep-ph/9703298 . 



R. J. Hill and M. Neubert, \Nucl. Phys. B657 (2003) 229-256 
|arXiv : h ep-ph /021 1018| 



C. M. Arnesen, J. Kundu, and I. W. Stewart, Phys. Rev. D7 2 (2005) 114002 
|arXiv : hep J ph7 05082 14} 

M. B. Voloshin, \Nucl. Phys. B154 (1979)~365J 



M. Beneke, A. Signer, and V. A. Smirnov, \Phys. Lett. B454 (1999) 137-146 
arXiv : hep-ph/9903260 . 

V. S. Fadin and V. A. Khoze, JETP Lett. 46 (1987) 525-529. 



A. H. Hoang et al, Eur. Phys. J. direct C2 (2000) 1, arXiv:hep-ph/0001286 



R Artoisenet, F. Maltoni, and T. Stelzer, JEEP 02 (2008 ) 102 
larXiv: 0712. 27701 I hep-phi I 



M. Czakon and A. Mitov, Phys. Lett. B680 (2009) 154-158, 
arXiv: 081 2. 03531 I hep-phi I 



[82; 
[83; 

[84] 
[85] 

[86] 
S71 



N. Brambilla, A. Pineda, J. Soto, and A. Vairo, Rev. Mod. Phys. 77 (2005) 1423 
|arXiv : hep-ph/ 04 10047] 

E. H. Wichmann and C.-H. Woo, 

Journal of Mathematical Physics 2 (1961) no. 2, 178-180 



http : //link . aip . org/link/? JMP/2/178/1 . 

A. H. Hoang, A. V. Manohar, I. W. Stewart, and T. Teubner, 
\Phys. Rev. D65 (2002) 0140141 |a7Xiv:hep-ph/0107144[ 



A. Pineda and A. Signer, \Nucl. Phys. B762 (200 7) 67-94J |arXiv : he p-ph/0607239 . 
W. Beenakker, R. Hopker, M. Spira, and P. M. Zerwas, 



Nucl. Phys. B492 (1997) 51-103, arXiv:hep-ph/9610490 



I. I. Y. Bigi, V. S. Fadin, and V. A. Khoze, Nucl. Phys. B377 (1992) 461-479 



W. Beenakker et al, arXiv: 1006.477 11 fhep-phl I 

A. D. Martin, W. J. Stirling, R. S. Thorne, and G. Watt, 



\Eur. Phys. J. C63 (20097T89- 285, a rXiv: 0901 .00021 I h ep-phi , 



[89] 
[90] 



P. M. Nadolsky et a/., Phys. Rev. D78 (2008) 013004, arXiv: 0802. 0007 [hep-ph] 



J. C. Collins and D. E. Soper, Nucl. Phys. B194 (1982) 445 



70 



[91] T. Becher and M. D. Schwartz, JEEP 02 (2010) 040, arXiv: 0911 .0681 [hep-ph] . 
[92] J. T. Goldman and H. Haber, Physica 15D (1985) 181-196. 



[93] Y. Kats and M. D. Schwartz, JEEP 04 (2010) 016, arXiv: 0912. 0526 [hep-ph] 



71 



